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ABSTRACT 

Galaxy clusters provide powerful laboratories for the study of galaxy evolution, particularly the 
origin of correlations of morphology and star formation rate (SFR) with density. We construct visible 
to MIR spectral energy distributions (SEDs) of cluster galaxies and use them to measure stellar masses 
and SFRs in eight low redshift clusters, which we examine as a function of environment. A partial 
correlation analysis indicates that SFR depends strongly on R/R200 (> 99.9% confidence) and is 
independent of projected local density at fixed radius. SFR also shows no residual dependence on 
stellar mass. We therefore conclude that interactions with the intra-cluster medium drive the evolution 
of SFRs in cluster galaxies. A merged sample of galaxies from the five most complete clusters shows 
(SFR) cx (i?/i?20o) 13±0 ' 7 for galaxies with R/R200 < 0.4. A decline in the fraction of SFGs toward 
the cluster center contributes most of this effect, but it is accompanied by a reduction in SFRs among 
star-forming galaxies (SFGs) near the cluster center. The increase in the fraction of SFGs toward 
larger R/R200 and the isolation of SFGs with reduced SFRs near the cluster center are consistent 
with ram pressure stripping as the mechanism to truncate star formation in galaxy clusters. We 
conclude that stripping drives the properties of SFGs over the range of radii we examine. We also 
find that galaxies near the cluster center are more massive than galaxies farther out in the cluster 
at ~ 3.5a, which suggests that cluster galaxies experience dynamical relaxation during the course of 
their evolution. 

Subject headings: galaxies:clusters:general, galaxies:evolution, galaxies:star formation, infrared radia- 
tion 
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1. INTRODUCTION 

The current paradigm for the evolution of the uni- 
verse and the growth of structure is based largely on 
observations of the luminous matter in the universe, i.e. 
individual galaxies, groups and clusters, and the cos- 
mic microwave background. While galaxy formation 
physics could previously be neglected, the era of precision 
cosmology increasingly demands detailed knowledge of 
galaxy formation to map obser vations of luminous mat - 
ter onto dark matter halos (e.g. Ivan Daalen et aT1l20TTI ). 
To do so precisely, we must understand the relationship 
between galaxy evolution and environment. 

Galaxy formation theory dates to the middle of the 
twentieth century. Early work explored t he physi- 
cal p rocesses responsible for star-formation (| Whipple! 
I1946D. speculated about the origins of the Milky Way 
( Egg en et al.l [19621 . and exami ned the impact of env i- 
ronment on galax y evolution (jSpitzer fe Baadd I1951D . 
lOsterbrockl ([1960D discovered that star-forming galax- 
ies (SFGs) are less common in galaxy clusters than 
in lower density environments, and this resu lt was 
subse quently re-exam ined with larger samples (|Gislerl 
119781: iDressler et all 119851 ). The dearth of vigorous 
star formation in galaxy clusters is mirrored by an 
under-abundance of spiral galaxies in these high den- 
sity regions, known as the morphology-density relation 
(|Dresslerl Il980t iPostman fc Gelled fl98l IDressler et all 
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fl997t [Postman et alJl200l . 

The impact of environment on the frequency and 
intensity of star-formation has been studied intensely 
in galaxy clusters and also at a variety of other density 
scales. These mea surements have employed both visible 
wave l ength colors dKodama fe Bower 2001; Balogh ~et al.l 
12001 IBarkhouse et al.l 120091: IHansen et al.l 120091) and 
emiss i on lines dAbraham et alJll996t IBalogh et al.lll997l 
20001 IKauffmann et al.l 12004: IChristlein fe Zabludofl 
2005t iPoggianti et alJ I200a iVerdugo et al l l200a 
Bragli aet al.l l2009t Ivon der Linden et al. | 20101) as well 
as mid-infrared (MIR) luminosities ( Bai et ahl 120061 : 
ISaintonge et~aT1l200a iBai et al.l 120091) . SFGs are con- 
sistently found to be more common and to have higher 
SFRs in lower density environments and at higher 
redshifts (jKauffmann et al.l 12004 IPoggianti etall 120061 
2008). This trend appears to reverse by z ~ 2, with 
star f ormation more c o mmon in clusters than in the 
field (jTran et al.l [2010l ; IHatch et al.l IMlTl) . However, 
even high-z cluster ga laxies form their sta rs earlier than 
coeval field galaxies ()Rettura et al.1 [20111 ). which is an 
expression of the so-called "downsizing" phenomenon 
HCowie et al.lU996h . 

The relationships between SFR, morphology and en- 
vironment in the local universe place strong constraints 
on models for galaxy evolution. Another important fac- 
tor is the presence of an evolutionary trend for galax- 
ies to have higher SFRs at higher redshifts. This was 
originally reported as an excess of blue cluster mem- 
bers at z ~ 0.4 compared t o z = (|Oemfed 119741 ; 
IButcher fc Oem ler 1978L IT984T) . and is commonly known 
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at the Butcher-Oemler Effect. This trend is now under- 
stood to track the simultaneous increase in the fraction 
of star-forming galaxies (SFGs) and in the SFRs of indi- 
vidual SFGs. An analogous trend has been examined 
in the MI R (iSaintonge et all 120081: lHaines et all 120091 : 
iTran et all 120101 : lHatch et al.ll20Tlf ). which is sensitive 
to dust-enshrouded star formation. 

The observed trends in star formation with en- 
vironment and the variation of these trends with 
redshift are usually attributed to changes in the 
sizes of cold gas reservoirs. Several mechanisms 
have been proposed to reduce galaxies' cold gas 
supplies and transform them from SFGs to passive 
galaxies. These mechanisms include ram-p r essure 
stripping of co l d gas (RPS: iGunn fc Gottl U972T: 
Abadi et al1IT99 9: Quili s et alJl200Ct IRoediger fe Hensler l 
2005] IRoe diger fe Bruggenl 120061 I2007t Ijachvm et al 
2007), ~gas starvati o n ([Larson et al.lll98d: IBalogh et all 
2000; iBekki et all 120021: iKawata fe Mulchaevi 120081: 



McCarthy et al.ll200 8: Book & Benson 20 101), galaxy ha - 
rassment ([Moore et all 119961 119981: lLake etaT] 119981), 
and i n teract i ons with the cluster tid al potential (jMerrittl 
119831 119841: INataraian et all 11 9981 ). Gas starvation 
operates throughout clusters, and it converts galax- 
ies from star forming to passive on a gas exhaustion 
timescale, which, f or normal spiral galaxies, is ~ 2.5 Gyr 
(Bigic i~e"t al.l l20lTT) . This time is similar to the cluster 
crossing time of 2.4 Gyr, which is the timescale appro- 
priate for dynamical processes like galaxy harassment. 
These timescales contrast sharply with the timescale ap- 
propriate for RPS, which truncates star formation on a 
gas stripping timescale, which is of order 10 5 yr. The effi- 
ciency of RPS also scales with ICM density, so it operates 
much more strongly near cluster centers than either star- 
vation or harassment. These differences mean that how 
rapidly star formation declines in cluster galaxies relative 
to the field constrains the mechanism primarily respon- 
sible for the removal of cold gas from cluster galaxies. 

The variation of SFR with environment can probe the 
relative importance of different environmental processes, 
but the conclusions drawn from apparentl y similar ob- 
servat ions sometimes conflict. For example. iMor an et all 
(|200l identified passive spirals in a sample of z ps 0.5 
clusters and determined that spiral galaxies rapidly turn 
passive when they enter the c luster env i ronme nt and 
then evolve into SO galaxies. iBai et all ()2009l ) argue 
that the similarity of the 24/xto luminosity functions 
observed in galaxy clusters and in the field suggests 
that the transition from star-formation to quiescence 
must be rapid, which implies that ram pressure strip- 
ping ( RPS) is the dominant mech a nism IVerdugo et all 
(2008) and Ivon der Linden et al.l (|2010f ). by contrast, 
find a significant trend of increasing SFR with radius 
to at least 2i? 2 oo from cluster centers. Because the trend 
of SFR with radius appears to extend to the virial shock 
(| White fc Frenldfl99lh . lv"on der Linden et al.l (|20Tot) con- 
clude that prepro cessing at the group scale is important. 
iPatel et all fj2009f ) find a similar trend for increasing av- 
erage SFR with decreasing local density down to group- 
scale densities (£ ga z ~ 1.0 Mpc~ 2 ). 

Evidence for pre-processing in groups is important be- 
caus e RPS is inefficient in low density gas so pr eprocess- 
ing (|Zabludoff fc Mulchaevi [1991: iFuiital l200l is likely 
driven by processes like gas starvation that operate in 



less dense environments. While preprocessi ng appears to 
be im portant in some groups and clusters, iBerrier et all 
(|2009l) found that very few cluster galaxies have previ- 
ously resided in groups, so the impact of preprocessing 
on a typical cluster galaxy must be limited. 

In Atlee et al. (2011; henceforth Paper I), we devel- 
oped a multi- wavelength method to identify AGNs based 
on either their X-ray luminosities or the shapes of their 
visible-MIR spectral energy distributions (SEDs). In this 
paper, we will employ this method to correct for AGNs 
and measure star formation in cluster galaxies. We use 
the results to study the relationship between star for- 
mation and the cluster environment. In particular, we 
consider the constraints placed on the important envi- 
ronmental processes that operate in clusters by the dis- 
tribution of star formation among cluster members. 

The paper is organized as follows: In $21 we review 
our observations, which are discussed in more detail in 
Paper I. In £[3] we review the mathematical formalism as- 
sociated with partial correlation analysis. In Sj4]we derive 
completeness corrections for the observed cluster mem- 
bers. We discuss the derivation of total infrared (TIR) 
luminosity functions (LFs) in jjH and in Sj6]we detail the 
results of our measurements. Finally, we examine the 
implications of these results for the environmental de- 
pendence of galaxy evolution in and we summarize 
these conclusions in |£8j Throughout this paper we adopt 
the WMAP 5-year cosmology — a A-CD M universe with 
n m = 0.26, n A = 0.74 and h = 0.72 (jDunklev et all 
l2009h . 

2. OBSERVATIONS AND MEMBER DESCRIPTION 

Paper I provides details of our photometry and spec- 
troscopy. It also develops methods to reliably identify 
low-luminosity AGNs and to measure galaxy properties 
like stellar mass and SFRs for identified cluster members, 
including AGN hosts. We briefly summarize the salient 
points below. 

2.1. Observations 

We identifie d cluster member gal axies using redshifts 
determined by iMartini et all (|2007|) . We supplemented 
these with redshifts from the literature, which we ob- 
tained from the NASA Extragalactic Database (NED) 2 . 
(These redshifts come from a variety of sources with un- 
known selection functions and success rates. See <JU) 

We have visible-wavelength images from the du Pont 
telescope at Las Campanas observatory, MIR images 
from the IRAC and MIPS instruments on the Spitzer 
Space Telescope, and X-ray images extracted from the 
Chandra archive. We measure visible, MIR and X- 
ray fluxes in redshift-dependent photometric apertures 
that approximate a fixed metric size. The aperture 
fluxes are then corrected to total fluxes at constant color 
with the -R-band Kron-li ke magnitude from SExtractor 
(jBertin fc Arno uts 1996). Our photometry spans the 
peak of the stellar continuum, so we have robust pho- 
tometric redshifts, which can identify catastrophic errors 
among the spectroscopic redshifts. We found 12 such 
catastrophic errors, one of which is an AGN (Paper I). 
We exclude these objects from our analysis. 

2 http:/ /nedwww. ipac.caltech.edu/ 
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2.2. Cluster Member Description and AGN 
Identification 

In Paper I, we constructed spectral energy distribu- 
tions (SEDs) from the photometry described in ^2.11 and 
we fi t models to t hese fluxes with SED template codes 
from lAssef et al.l (|2010D . These model SEDs are used 
to derive photometric redshifts and if-corrections. We 
also employed the models to measure MIR color correc- 
tions, and we used the results to determine rest-frame 
luminosities. 

We employ rest-frame, visible wavelength colors to de- 
termine the mass to light ra tio of each galaxy in the 
sample (|Bell k, de Jong! 120011 ). and we combine these 
with the measured luminosities to infer A/*. To mea- 
sure SFRs, we empl oy both 8pm and 24/im luminosi- 
ties fZh u et alJ[20(M ). When SFRs can be measured in 
both bands, we take the geometric mean of the two mea- 
surements. We found in Paper I that SFRs determined 
independently from 8 pm and 2Afim luminosities show a 
scatter of ~ 0.2 dex with respect to one another. This 
scatter reflects the systematic uncertainty in SFR mea- 
surements determined from either band separately, and 
the implied uncertainty is comparable to the systematic 
uncertainties in the measured stellar masses ( 0.3 dex, 
Paper I). 

Before we measure stellar masses and SFRs of cluster 
members, we identify and correct for AGNs. We employ 
two independenty methods to identify AGNs: the shapes 
of the model SEDs (IR AGNs) and the X-ray luminosi- 
ties measured with Chandra (X-ray AGNs). The hosts 
of IR AGNs are corrected for the presence of the AGN 
before we calculate M* and SFR. X-ray AGNs without 
visible signatures in their SEDs do not contribute signif- 
icantly to the measured visible and MIR fluxes, so we do 
not correct those objects. In Paper I, we found that the 
IR and X-ray AGN samples are largely disjoint. This im- 
plies that X-ray only AGN selection can overlook a large 
fraction (~ 35%) of AGNs. We explore the potential 
consequences of this bias in §6.51 and S 36.61 

3. PARTIAL CORRELATION ANALYSIS 

When confronted with a system of mutually correlated 
observables, it can be difficult to establish which vari- 
ables drive the correlations. Partial correlation anal- 
ysis measures the relationship between two variables 
with all other parameters held fixed and can identify 
which variable(s) control the observed correlations. Par- 
tial correlation analysis has been applied in the past 
to develop a fundam ental plane of black hole activity 
(jMerloni et al.l 120031 ) and to probe the dependence of 
SFR on both stellar mass and en vironment simultane- 
ously (IChristlein fc Zablud off 2005). We will use the sim- 
plest formulation of partial correlation analysis, which 
relies only on direct measurements and does not account 
for upper limits. 

The simplest case is a system of only three variables, 
Xj. This is called the first-order partial correlation prob- 
lem. The correlation coefficient for Xi and X2 at fixed X3 
can be expressed as 



cient (e.g. the Pearson or Spearm an coefficients) between 
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where pij is the standard two- variable correlation coeffi- 
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Higher order problems 



describe systems with more variables. For a system of N 
variables, the (N — 2) th order partial correlation coeffi- 
cient ry.i...jv\{y} of variables X4 and Xj can be written, 
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where dj = {-l) l+j M^ (|Kendell fc StuartlH977l) . M (J 
is a reduced determinant of the correlation matrix R, 
where Rij = p^, and is the two- variable correlation 
coefficient of Xi and Xj. The determinant Mij can be 
interpreted as the total correlation among the variables 
of the system in the absence of i and j. It is calcu- 
lated from R with the i t h row and j th column eliminated 
([Kendell fc Stuarti [19771) . 

Given a partial correlation coefficient from Eq. [2 we 
would like to know its significance. This can be evaluated 
from <Tij.i,..N\{ij}, 



1 



Hj.l...N\{ij} 



rij.i...N\{ij} 



y/m-N 



(3) 



where is the partial correlation coefficient 

given by Eq. [21 N is the number of variables in the sys- 
tem, and m is the number of objects in the sample. The 
statistical significance of Tij,i...N\Uj} i s determined from 
the Student's t-distribution with dispersion 0-^1 ...N\Uj} 
(|Wall fe Jenkinsl[200ll) . 

Partial correlation analyses can take both parametric 
and non-parametric forms. These are analogous to the 
more commonly applied two-variable correlation analy- 
ses. Equation [5] can be applied to any of the correlation 
coefficients in common use. However, Eq. [3] is defined 
for the parametric Pearson's correlation coefficient, so it 
is appropriate only for that estimator or the closely re- 
lated, non-parametric Spearman coefficient. We want a 
non-parametric approach, so we rely on Spearman cor- 
relation coefficients in our analysis. 

4. COMPLETENESS CORRECTIONS 

We wish to examine the distributions of the stellar 
masses and SFRs described in £12.21 to probe the impact 
of the cluster environment on galaxy growth. However, 
to do this we must first correct for selection effects. The 
spectroscopic selection function that defines our sample 
is unknown, because many of the sources that contribute 
to the redshifts in the literature do not define their tar- 
get selection functions or rates of success. Furthermore, 
the MIR observations do not uniformly cover the clus- 
ter fields. Therefore, we empirically determine both our 
spectroscopic and MIR selection functions to correct for 
these effects. 

4.1. Spectroscopic Completeness 

We examine only spectroscopically-confirmed cluster 
members. Many o f the redshifts we use come from 
iMartini et al.l (|2007l ). which we supplemented with red- 
shifts from other sources in the literature. This results in 
a complex selection function that is poorly known a pri- 
ori. However, this completeness function is required to 
correct the properties of observed cluster galaxies to the 
intrinsic distribution for all cluster members. We take 
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an empirical approach to determine spectroscopic com- 
pleteness and correct the measured cluster members to 
the total cluster galaxy population. 

For each cluster, we bin galaxies identified in the pho- 
tometric source catalog by V — R color, i?-band magni- 
tude and i?/i?20o- We hnd significant variations in the 
fraction of galaxies with spectra (f spe c) as a function of 
i?/i?200 and tur, but the variation with color is at most 
minor. A partial correlation analysis of f spe c as a func- 
tion of color, magnitude and position shows no signifi- 
cant partial correlation with V — R at 95% confidence 
in any cluster, while / spec correlates with both m,R and 
R/R200 at > 99.9% confidence. We therefore collapse 
the measurement along the color axis and determine the 
fraction of galaxies with spectroscopy as a function of 
-R-magnitude and position only. This results in better 
measurements due to the larger number of galaxies per 
bin. 

The fspec measured above is one way to express 
the spectroscopic completeness of galaxies in a given 
magnitude-radius bin. However, what we really want is 
an expression for the spectroscopic completeness, C spec , 
of cluster members, 



r _ N C l,spec( x) 

^ spec \?b ) 



N C i(x) 

N spec 

Ntot{x) Nspecix) ~ N C i(x) 



Net 



,spec(x) N tot (x) 



(4) 
(5) 



where x is the position of a given bin in magnitude-radius 
space, Nci, S pec is the number of galaxies with spectra 
that are cluster members, Nci is the number of true clus- 
ter members, N spe c is the number of galaxies with spectra 
in the cluster field, and N to t is the number of galaxies in 
the input catalog. All of these quantities except Nci can 
be measured directly from the input catalogs. We would 
need to infer Nci using some additional piece of infor- 
mation, so we prefer to rely on f spe c rather than C spec if 
possible. 

If the redshifts reported in the literature were not pre- 
selected for cluster membership or if the redshift fail- 
ure rate was high, f spec (x) would be a good proxy for 
C sp ec(x), and the approach in Eq. 2] would be unneces- 
sary. If this were the case, the fraction of galaxies with 
spectra that are cluster members (f mem ) should drop 
with R/R200 as the fraction of field galaxies increases. 
Figure Q] shows that / mem does not always trace the de- 
cline in the density of cluster galaxies. This implies that 
fspec is not a good tracer of C spec , and the more sophis- 
ticated approach of Eq. 0]is required. 

Before we can employ Eq. [4j we need to know the num- 
ber of cluster galaxies in each bin. To do this, we esti- 
mate the number of field galaxies in the bin with the 
i?-band magnitude-numbe r density relation reported by 
iKummel fc Wagnerl (|2001f ). We subtract the field galax- 
ies from the total number of galaxies in the bin to esti- 
mate the number of cluster galaxies. 

This approach introduces two types of uncertainty. 
The first is simple Poisson counting uncertainty due to 
the small number of field galaxies, typically a few to 10 , 
in each bin. The second is cosmic variance. lElhsl (fl987h 
reports a £>-band magnitude-number relation that in- 
cludes measurements from a number of other authors. 
The different surveys use fields of different sizes, so the 
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Fig. 1. — Comparison of galaxy density (upper panel) and spec- 
troscopic membership fraction (lower panel) as a function of radius 
for the clusters in our sample with enough confirmed members to 
make a meaningful measurement. If / spec were a good proxy for 
fci,speci the upper and lower panels would have similar slopes. 
None of the clusters in the sample exhibit this behavior. 

scatter of their results about the best-fit relation provide 
a measure of the cosmic variance, which contributes of 
order 10% uncertainty on the number of field galaxies in 
a typical bin. The number of field galaxies in a given bin 
depends on magnitude and cluster mass, but it gener- 
ally ranges from 1-10 galaxies. At faint magnitudes, the 
number of field galaxies is generally comparable to the 
number of cluster galaxies, and Poisson fluctuations in 
the number of field galaxies dominate the uncertainties 
in the completeness measurements. The spectroscopic 
completeness measurements and associated uncertainties 
for each cluster are summarized in Table [TJ 

Figure [2] shows the spectroscopic completeness (C spec ) 
for 6 of the 8 galaxy clusters in our sample. The re- 
maining 2 clusters (A644 and A2163) have too few con- 
firmed members to make a reliable measurement. The 
dashed, vertical lines on the right column of Figure [2] 
indicate the observed magnitude that corresponds to 
Mr — — 20 for the average K-correction in each cluster. 
T he follow-up spectros copy of X-ray sources conducted 
bv lMartini et al.l f|2006T ) is only complete to this luminos- 
ity limit. Clearly, completeness becomes quite poor for 
Mr > —20 in all clusters, so we restrict our sample to 
galaxies with Mr < — 20. 

We also considered completeness as a function of lumi- 
nosity and stellar mass instead of itir. However, these 
quantities have higher uncertainties than observed mag- 
nitudes, especially for galaxies without spectroscopic red- 
shifts to fix their distances. Therefore, we measure com- 
pletness as a function of tur and i?/i?2oo- 

4.2. Mid-Infrared Completeness 

The depth of the MIR images varies as a function of 
position across the clusters. This is a result of the Spitzer 
mosaicking schemes, which were chosen to provide good 
coverage of the known X-ray point sources in the cluster. 
These mosaic schemes lead to variations in the number 
of overlapping images, and therefore to variations in sen- 
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Fig. 2. — Fraction of cluster galaxies with spectra as a function 
of projected distance from the cluster center (left) and i?-band 
magnitude (right). The measurements in each column have been 
separated by the indicated variable (colored points). The black 
points show the completeness averaged over all of the colored bins, 
and the grey bands show the la confidence intervals on the total 
completeness in each cluster. The dashed, vertical lines on the 
righthand column indicate the magnitude corresponding to Mr 
—20 mag. The average measurements use only galaxies with Mr < 
—20 mag and R/R200 < 0.4 (filled, black points). 

sitivity, across the cluster fields. 

In addition to these sensitivity variations, the Spitzer 
footprint features some non-overlapping coverage by the 
different IRAC bands. This results from the IRAC map- 
ping strategy, which simultaneously images two adjacent 
fields in different bands. The pointings chosen by the ob- 
server then determine the degree of overlap between the 
IRAC channels. For a galaxy to enter the final sample, 
it must include detections in at least 5 bands to ensure 
that the fit results for that galaxy are well constrained. 
This means that a faint galaxy in a region of a cluster 
with overlapping 3.6/im and 4.5/xm images, for example, 
might be more likely to appear in the final sample than 
an identical galaxy in a part of the cluster with only 
4.5/ito coverage. 

To construct ensemble statistics for whole clusters, 
we require sensitivity corrections that account for vari- 
able depth across the cluster fields and for the differ- 
ent footprints in the Spitzer bands. We again take 
an empirical approach to completeness correction. We 
measure the MIR flux uncertainties at the locations of 
all confirmed cluster members from the Spitzer uncer- 
tainty mosaics. A t each position, we combine the two 
lAssef et al.1 (|2010T l star-forming templates with arbitrary 
flux normalizations 1000 times to produce galaxies with 
10~ 2 < SFR/1 M Q yr- 1 < 10 2 . From these artificial 
galaxy SEDs, we construct model fluxes and determine 
whether the galaxy represented by each SED would have 
been detected at 3<r based on the flux uncertainty at 
each position. We bin the results by flux and by i?/i?200 
to estimate completeness separately at 8/im and 24/im. 
Figure [3] shows the results of this measurement for the 
6 clusters in Figure [5] IRAC and MIPS completenesses 
clearly depend on both flux and i?/i?2oo- 

The uncertainties in MIR completeness result from the 
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Fig. 3. — MIR completeness as a function of flux for the 8fim (left) 
and the 2Afim (right) Spitzer bands. In each column, the sample 
has been separated into 4 radial bins. Fluxes have not been color- 
corrected and are given in the observer frame. Uncertainties are 
shown for a single radial bin to indicate typical values. MS 1008.1- 
1224 was not observed wit h M IPS. Completeness measurements 
are derived as described in 34.21 

incomplete spectroscopic sampling of the galaxies in a 
given bin. We implicitly assume that the identified clus- 
ter members in each bin are representative of the be- 
havior of the unidentified members. This assumption 
means that the precision of the completeness correction 
in a given bin is fixed by the number of identified clus- 
ter members in that bin. The completeness we measure, 
Cmir, is our best estimate of the "true" MIR complete- 
ness C M U IR associated with a hypothetical, spectroscopi- 
cally complete sample of galaxies. In a bin with N clus- 
ter members, the expected number of detections is sim- 
ply C M ™ R N. However, the actual number of detections 
will have some range around Cj^"j R N, which leads to an 
uncertainty in the inversion of Cmir to a completeness 
correction. This uncertainty is set by the expected vari- 
ation in the number of galaxies, which is best described 
by binomial statistics. This allows calculation of asym- 
metric error bars on Cmir and accounts naturally for 
upper and lower limits. Typical uncertainties returned 
by this procedure are ~ 20%. The full set of MIR com- 
pleteness measurements and the associated uncertainties 
are summarized in Table [5J 

4.3. Merged Cluster Sample 

We have defined the MIR completeness measurements 
in Figure [3] so they apply only to galaxies with spectro- 
scopic redshifts. The two corrections, applied serially, 
give total completeness corrections. The total correction 
Xq for a galaxy G is, 



X G = 



1 



1 



C spec {RG/ RwQ,mR,G) 



Cmir(Rg/R2oo, Sv,g) 
(6) 



where C spec is the spectroscopic completeness (Figured]) 
and Cmir is the MIR completeness (Figure[3]) . The com- 
pleteness corrections described by Eq. [5] can be applied 
to individual galaxies to extrapolate from the measured 
galaxy samples to the full cluster population. In cases 
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where multiple corrections can be derived for a single ob- 
ject, we combine these corrections in the same way the 
data are combined. For example, the completeness cor- 
rection for a galaxy with SFR measurements from both 

8fim and 24/im fluxes is given by X — [Xs fml X24 l mi] 1 ^ 

because SFR = [SFR^ m SFR 2i , m ] 1/2 (gpj). 

To examine the dependence of star formation and black 
hole growth on environment, we need to construct a 
merged cluster galaxy sample. We identify 5 clusters 
(A3128, A2104, A1689, MS1008 and AC114) with the 
best completeness estimates and combine their members. 
The relatively small number of galaxies in A3125 results 
in highly irregular behavior of the completeness func- 
tions. As a result, any corrections we applied would 
depend critically on the binning scheme. We therefore 
exclude it from the main cluster sample. 

The clusters in the main sample can be stacked to yield 
better statistics. To construct the stacked cluster, we 
weight individual galaxies by their completeness correc- 
tions (Xg)- The correction is a combination of the spec- 
troscopic and photometric completeness corrections from 
gZDand S3 as given by Eqs. HE 

5. LUMINOSITY FUNCTIONS 

Luminosity functions (LFs) provide an important di- 
agnostic for the difference between cluster galaxies and 
field populations, because LFs are sensitive to the en- 
tire clust er populati o n rath er than only the average. For 
example, iBai et aLl (|2009( ) employed the total infrared 
(TIR) LF to infer that RPS controls the evolution of 
SFRs in cluster galaxies. In this section, we discuss the 
derivation of total infrared (TIR) luminosities and our 
method to construct luminosity functions. We discuss 
the results in §6.41 

5.1. Total Infrared Luminosity 

Our MIR observations cover a relatively narrow wave- 
length range from 3.6fim to 24/ito. To compare our re- 
sults with previous studies, we need to infer Ltir from 
the observed Ls^ m and L24 fim , so we must apply bolo- 
metric corrections (BCs). To esti mate Lttr from the 
Spitzer luminosities, we employ the lDale k. Heloul ()2002f ) 
SED template library, which includes a wide variety of 
SEDs. These SEDs differ from one another according 
to the parameter a, which describes the intensity of the 
radiation field on a typical dust grain. 

Before we calculate BCs for IR AGNs, we first sub- 
tract the AGN cont ribution (Paper I). We then fit each 
iDale fc Heloul ()2002D template to the rest- frame 5.8, 8.0 
and 24/ito fluxes and use the template that best fits the 
data to measure our fiducial BCs. In the frequent cases 
where luminosities in one or more of these bands are 
unavailable, we estimate the m issing luminosities from 
model SEDs (|Assef et al.ll2010h . When this is necessary, 
we assign uncertainties to the model fluxes from the un- 
certainties on the model SED. We only calculate Ltir 
for galaxies with detections in at least one of the 8^m 
and 24/iTO bands. 

In galaxies that have measurements of both Lg^m and 
we calculate Ltir separately for each band and 
take the geometric mean of the results. This follows our 
treatment of SFRs in Paper I. In other cases, we sim- 
ply use the BC appropriate for the band where we have 



a detection. Typical BCs are ~ 6 for Lsum and ~ 8 
for 

L24fim- We also construct 68% confidence intervals 
for each BC based on the Ay 2 = 1 interval for each 
galaxy. These uncertainties are asymmetric, and they 
add in quadrature to the uncertainties on Lg and L24 to 
give the total uncertainty on Ltir- 

5.2. Luminosity Function Construction 

We construct luminosity functions from galaxies in 
the main cluster sample whose luminosities — Ltir. — 
L~^t — we determine as described in £15.11 If we account 
for the uncertainties on Ltir, we can reduce our sensi- 
tivity to Poisson fluctuations in the number of luminous 
galaxies. We distribute the galaxy weights described 
in §3 over luminosity bins according to the probability 
that the true luminosity of a galaxy with best-estimate 
L = Ltir lies in a given bin. Due to the uncertainty 
on the LF prior, this technique increases the statistical 
uncertainty on the total weight in each bin by ~ 10%. In 
exchange, we reduce the much larger uncertainty intro- 
duced by stochasticity in the number luminous galaxies. 

To distribute galaxy weights over luminosity bins, 
we employ an asymmetric probability density function 
(PDF) that considers o\ and a u separately. We inte- 
grate the PDF across each luminosity bin to determine 
the weight in each bin, which we add to construct the 
total LF. The PDF we employ is piecewise smooth, and 
it approaches the Normal distribution when a u ~ o - ;. It 
is described in more detail in Appendix [A} 

In addition to the PDF, we require a prior on the shape 
of the LF to correct for Eddington-like bias due to the 
steepness of the LF above L*. We a dopt a Schechter 
function fit to the Coma cluster LF from lBai et al.l (2006) 
as the baseline prior. We then correct the Coma LF to 
the redshift of each individual cl uster according to th e 
evolution of the field galaxy LF (|Le Floc'h et alJ feOOS). 
We add the uncertainty on the prior to the statistical 
uncertainty on the LF in each luminosity bin. The prior 
has a strong impact on the bright-end shape of the LF 
because there are few cluster galaxies to constrain the 
LF in this regime. The results are discussed in §6.41 

6. RESULTS 

We apply the weights derived from the completeness 
corrections described in §J] to the main cluster sample, 
which is a subset of the cluster galaxies in Tabled These 
corrections allow us to examine the average environmen- 
tal dependence of M* (§M|) and SFR f §Ol6~5)) . and the 
redshift dependence of star formation f H6.6|) . Before con- 
ducting these analyses, we perform a partial correlation 
analysis to determine which observed properties of galax- 
ies in clusters most strongly correlate with star formation 
( §6.ip . The results inform the rest of our work. 

6.1. Partial Correlation Analysis 

The cluster environment can significantly alter the 
evolution of cluster member galaxies, as described in 
§T] However, when we attempt to distinguish be- 
tween the physical properties that might cause these ef- 
fects, we confront a system of mutually-correlated ob- 
servables. For example, SFR depends on both pro- 
jected lo cal galaxy density to the 10 th nearest neigh- 
bor (E in : I0sterbr^[l96a [Qemlerj[l97l lDressledll980t 
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iKauffmann et all I2004D and position within the clus 
ter (i?/-R 2 on:lKodama fc Bowerll200ll;lBalogh et alll2004 ; 



Christlein fc Z abludoff 120051 iBlanton fc Berlindl 12007 
Hansen et al J 1200a Ivon der Linden et al.l l2010f) . Figure 



H demonstrates the correlations between SFR, position, 
projected galaxy density and M* among SFGs. It is not 
immediately clear which of these is the most fundamen- 
tal. 

While the causal connection between morphology and 
the local density of galaxies is well established (e.g. 



Dressier 1980; Dres sier et alJll99alPostman et al.ll2005D 
Moran et al.l ([20070 find strong evidence that the mor- 



phologies and star-formation rates of massive spiral 
galaxies in clusters evolve separately. This implies that 
a factor other than local density may control star for- 
mation in cluster galaxies. Because M*, R/R200 and 
projected local density (£10) are all mutually correlated, 
it is not easy to determine which variable(s) drive the 
environmental dependence of star formation. Therefore, 
we use a partial correlation analysis to disentangle these 
dependencies. The mathematical formalism for partial 
correlation analysis is described in We do not con- 
sider completeness corrections for this analysis, so we 
include galaxies from all 8 clusters. 

We consider only objects with measurements of all pa- 
rameters under consideration and ignore galaxies with 
upper limi ts. This differs from th e simil ar analysis con- 
ducted by IChristlein k Zabludoffl (|2005l ). who also con- 
sidered upper limits. As a result, ou r resul ts are more 
sensitive than Christlein k ZabludofJ (|2005[) to system- 
atic effects like variations in sensitivity within or between 
clusters. Because of this, we do not rely directly on the 
strength of any partial correlations, but only on the pres- 
ence or absence of such correlations. For variables with 
significant partial correlations, we perform stacking anal- 
yses, which can account for incompleteness. (See ^6.21 



We perform a partial correlation study on a sys- 
tem of five variables: SFR, M*, R/R200, the 
iDressler fc Shectmanl (|1988[ ) substructure parameter (5) , 
and projected local density of cluster members (Sio)- 
The partial correlation coefficients returned by the anal- 
ysis are listed in Table 0] 

Table 2] shows that SFR depends strongly on i?/i?200 
(fs, partial = +0.34), but it shows no significant de- 
pendence of SFR on M* once the influence of R/R200 
has been factored out (rs, pa rtial = +0.09). This con- 
flicts with earlier results, whi ch generally find either 
that SFR depends o nly on M* (jGrutzbauch et al.|[20TTI: 
iRettura et all 120111) or that SFR depends on b oth M* 
and environment ( Christlein fc Zabludoff 2005). One 
reason for this discrepancy is that Table |4] does not in- 
clude non-SFGs. The fraction of non-SFGs is higher 
among more massive galaxies, which introduces a depen- 
dence of (SFR) on M*. Another possible factor is the 
difference in stella r mass es we examine. The sampl es of 
iGriitzbauch et al.l (|201lD and IRettura et all (|2011f l ex- 
tend down only to ~ 3 x 10 10 M Q , while our sample 
extends to approximately 10 10 M Q . If the cluster envi- 
ronment affects lower mass galaxies more strongly, we 
should be more sensitive to this effect. Because the effi- 
ciency of RPS should depend on M* , our expanded mass 
range improves our ability to distinguish between RPS 
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Fig. 4. — Correlations of star formation with position in the clus- 
ter (top row), projected local density (middle row), and stellar mass 
(bottom row). Galaxies with no measurable star formation are ne- 
glected. Colors denote the different clusters in the sample: A3128 
(black), A3f25 (red), A644 (blue), A2104 (green), A1689 (cyan), 
A2f63 (magenta), MS1008 (orange), AC114 (violet). Large black 
points show the median values of the galaxy sample after it has 
been binned by SFR. SFR shows strong correlations with both 
R/R200 an d S10 (rg = +0.35 and rg = —0.34, respectively), but 
no correlation with M*. Partial correlation coefficients derived 
from these data are listed in Table |4"1 



and gas starvation. 

Interestingly, SFR also shows no residual dependence 
on S10 (rs, P artiai = +0.02), which suggests that the 
SFRs of cluster members are driven by the local condi- 
tions of the ICM rather than by interactions with nearby 
galaxies. We emphasize that this conclusion applies to 
SFR only, and we do not consider the dependence of mor- 
phology on Sio or R/R^oq. This r esult is also consistent 
with results from previou s authors (|Poggianti et al.|[l999t 
IMoran et al.l [2006. 2007]) who found that the processes 
that alter SFR and morphology are likely to be physically 
distinct. The determination that SFR is more closely re- 
lated to R/R 200 than to £10 distinguishe s our results 
from those of IChristlein fc Zablu doff (|2005[) . who do not 
discriminate between different environmental tracers. As 
a result, they are agnostic about the process(es) that 
drive the (iS-P-fi^-radius relation. Our analysis also relies 
on SFGs alone, which disting uishes it from the work of 
IChristlein fc Zabludoffl (|2005t ). who included upper lim- 
its for galaxies with no measurable star formation. This 
accounts for the lack of a strong anti-correlation between 
M* and SFR, which is driven by a decline in the frac- 
tion of SFGs at higher M* rather than a reduction in the 
SFRs of individual SFGs. 

SFR also shows no relationship to local substructure, 
as measured by S, at fixed R/R200. Indeed, even a 
two-variable correlation test returns no correlation be- 
tween <5 and SFR (rg = -0 .03). This conflicts with 
IChristlein fc Zabludofl (|2005D . who reported a strong 
correlation of SFR with local substructure. However, 
S requires a robust spectroscopic sample from which to 
measure local velocity dispersions. As a result, the sub- 
structure measurements for some of our clusters with less 
complete spectroscopy are probably unreliable. We re- 
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Fig. 5. — Radial distributions of all cluster members scaled to 
R200- Panel (a) compares the radial distributions of low-mass 
(blue dashed) and high-mass (red solid) galaxies, divided into two 
equally-sized subsamples at M C ut = 3.9 X 10 10 Mq. The two dis- 
tributions differ at 99.9% confidence a fter we exclude BCGs as 
defined by von dcr Linden ct al. (2007). Panel (b) compares the 
radial distributions of galaxies with (blue dotted) and without (red 
dashed) an 8fim flux excess to the distribution of all galaxies with 
8/im detections (heavy black). The distribution of galaxies with 
no measurable excess shows a marginal difference compared to the 
distribution of all cluster members (95% confidence). 

peat the test with A3128 and A3125, which have the 
most complete spectroscopy and are the only two clus- 
ters with significant substructure. The results are indis- 
tinguishable from the full cluster sample. 

6.2. Mass-Radius Relation 

In M6. 11 we reported a strong correlation of SFR with 
R/R200 but no residual dependence of SFR on M* or M* 
on radius. This mi ght indicate that M» is inde pendent 
of environment, as Ivon der Linden et al.l (|2010T l found. 
However, it might also mean that the SFR-i?/i?200 cor- 
relation is strong enough to eclipse any more subtle cor- 
relations that might appear among a sample composed 
entirely of SFGs. The galaxy sample examined in £16.11 
includes only a few hundred galaxies, and the sample 
preferentially excludes the most massive galaxies, which 
tend not to show active star formation. As a result, 
HQ] might show no correlation between M* and R/R20Q, 
even if the full cluster gal axy sample includes one. 
iChristlein fc ZabludofJ (|2005l) found a strong partial cor- 
relation of mass with R/R20Q. This correlation would be 
difficult to produce if BCGs alone produce a false correla - 
tion of M* with R/R200, as Ivon der Linden et al.l (|2010f ) 
claim, because normal cluster galaxies are much more 
numerous than BCGs. 

To test whether our data support the presence of a 
radial trend in M*, we look for direct variations of M* 
with R/R200 without regard to correlations with other 
variables. We first divide the galaxy sample into two 
samples with equal numbers of galaxies, and we apply 
a K-S test to check for a difference between their ra- 
dial distributions. For this analysis, we include members 
of all 8 clusters, and we e xclude BCGs as defined by 
Ivon der Linden et al.l (|2007l ) from the sample. The re- 
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Fig. 6. — Average stellar mass for galaxies in the stacked cluster 
sample as a function of radius. Brightest cluster galaxies are ex- 
cluded from the fit because of their unusually large stellar masses 
and SFRs compared to other galaxies near the centers of clusters. 
The heavy line indicates the best-fit power-law to the 6 outer bins. 
The two innermost radial bins were excluded from the fit based 
on their large excesses. This resulted in a reduction in total x 2 
from 9.1 to 1.0. The best fit yields M* oc [R/R200] 0,4±a2 , and 
the shaded region indicates the 68% confidence interval to the fit. 
The residuals are shown in the lower panel. 

suits are shown in Figure [5^. The K-S test returns a 
probability < 0.1% that the high- and low-mass samples 
have the same radial distributions, so massive galaxies 
are preferentially found closer to the centers of their par- 
ent clusters, even in the absence of BCGs. We weight 
members of the main cluster sample by their complete- 
ness to determine the average mass as a function of ra- 
dius. The average mass in a given bin is, 

where N is the number of galaxies in the bin with Mr < 
—20. The Wi are the weights derived from Eq. H) Figure 
[5] shows the resulting average masses as a function of 
radius. 

The innermost radial bin Figure [6] shows a strong ex- 
cess compared to the other bins, and the second bin hints 
at an excess. We fit a power law to the six outer ra- 
dial bins in Figure [5] (solid line) to measure the strength 
of the mass excess in these bins and to determin e why 
our results differ from Ivon der Linden etaTI (|2010l ). The 
slight increase of (M*) with i?/i? 2 oo beyond R « 0.1i?2oo 
is consistent with the tendency of massive galaxies to be 
accreted more recently than less massive galaxies and for 
recently accreted gala xies to lie further out in the cluster 
(|De Lucia et aLlfeOllft . Alternatively, galaxies presently 
near the center of the cluster spend more time near the 
cluster center, on average, than galaxies farther away. 
Therefore, they are subject to stronger tidal forces from 
the cluster potentia l and will lose more o f their mass 
(|MerrittlH98l 11984 iNataraian et all 11998ft . The large 
uncertainties on the fit preclude any attempt to distin- 
guish between these scenarios. 

Figure [6] shows that the average masses of galaxies near 



Environmental Effects on SF 



9 



the cluster center show an excess compared to their coun- 
terparts further out in the cluster. The innermost ra- 
dial bin in Figure [6] differs from the best fit model by 
2.8<7, and the mass excess in the second radial bin is 
significant at 1.8a. This indicates that the cluster core 
(R < 0.05i?2oo) tends to host more massive galaxies than 
the outer regions, even if we neglect BCGs. Mass segre- 
gation among cluster galaxies can be introduced as the 
cluster relaxes to virial equilibrium. The dynamical re- 
laxation timescale in the inner mass bin for a cluster at 
z = 0.15 with a = 1200 km s _1 is approximately 300 
Myr, so the core of a typical cluster in our sample should 
be dynamically relaxed. For the same "typical" cluster, 
the crossing time for the sphere defined by R < 0.4i? 2 oo 
is approximately 900 Myr. This crossing time implies a 
dynamical relaxation time of 7 Gyr, which is longer than 
the age of a cluster at z = 0.15 that "assembled" at z = 1 
(5.8 Gyr), so the sphere with R < 0.4i?2oo has not yet 
had time to relax. 

Fi nally, FigureElsuggests a reas on for our disagreement 
with [van der Linden et all (|2010l) . The signal comes pri- 
marily inside 0.05j?2nn , which corresponds to ~ If 5 at the 
median redshift of the Ivon der Linden et al.l (j2010T ) sam- 
ple (z ps 0.8). Due to SDSS fiber collisions, only a few 
galaxies inside 0.05i?2oo will have redshifts in each clus- 
ter. This shifts the media n of th e innermost radial bin in 
the Ivon der Linden et al.l (|2010T ) sample to ~ 0.08i?2oo- 
This is comparable to the second radial bin in Figure [6l 
If this was our innermost bin, we would not find any de- 
pendence of M* on R/R^nn, so the disagreem ent between 
our results and Ivon der Linden et al.l (|2010l) likely result 
from fiber collisions in SDSS. 



6.3. Environmental Dependence of SFR 

In Paper I, we examined the R/R200 distributions of 
AGNs and found no significant difference between the 
positions of AGNs and normal cluster members. The 
lack of radial dependence among AGNs could due to the 
small sample size, it could indicate a weak dependence 
of the amount of cold gas on i?/i?200i or h might mean 
that AGN fueling is poorly correlated with the total cold 
gas reservoir of its host galaxy. To test these hypothe- 
ses, we define a sample of galaxies with 8/im flux excesses 
as those galaxies whose measured 8/im flux exceeds the 
flux expected from a passively evolving galaxy matched 
in Mj< at more than 2a. Figure [SJj compares the radial 
distributions of galaxies with and without an 8/im excess. 
These objects include both SFGs and AGNs. We again 
excluded BCGs from these samples. The radial distri- 
bution of galaxies with 8/im excesses is indistinguishable 
from the merged sample, but galaxies without an excess 
are located closer to the centers of their host clusters 
than the average cluster galaxy at 95% confidence. 

The dependence of dust emission on R/R200 
shown in Figure \5jp is consistent with the es- 
tablished depen dence of SFR on position within 
galaxy clusters dKodama fc Bowerl 120011: iBaloeh et al.l 



20041: iChristlein fc Zabludofl 120051: lHansen et all 120091 : 
von der Linden et al.l 20101 ) and with our results in §6.11 
One way to test the origin of this effect, and by extension 
the SFR-density and SFR-radius relations, is to measure 
the average SFR as a function of radius. We weight indi- 
vidual SFGs by their total completeness (Eq. [6J and bin 
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Fig. 7. — Average star formation in the stacked cluster sample 
as a function of position. Each panel shows two different binning 
schemes: fine (black squares) and coarse (orange triangles) for the 
same galaxy samples. The top panel shows average SFR among 
all galaxies. The solid line indicates the best-fit power law to the 

data (SFR oc [R/ifeoo] 1 3 ° 7 ), and the red line shows the gas 
starvation model of Book & Benson (2010), normalized to match 
the observed SFRs. The shaded region indicates the la confidence 
interval to the fit. The middle panel shows the fraction of SFGs 
(SFR > 3 Mq yr~ x ) as a function of position, with the best fit 
(fsF oc [R/-??20o] 0,4±0 ' 2 ) shown by the line and the la confidence 
interval shown by the shaded region. The bottom panel shows the 
averaged SFR among SFGs ((SFR)sfg) versus R/R200- Galaxies 
with R/R200 $ 0.1 have lower (SFR)gpQ than galaxies outside 
O.IR200 at > 99.9% confidence. 

them in radius to determine (SFR): 



(SFR) 



\WSF, 



iSFR] 



(8) 



where wsfa and Wj are the weights for SFGs and all 
galaxies, respectively. For this calculation, we define 
SFGs as all galaxies with SFR > 3 M Q yr' 1 . This guar- 
antees that we are not subject to biases due to variable 
sensitivities across the cluster fields. 

We fit a power law to (SFR) as a function of R/R200 
and find, 

log w [SFR] = (1.3±0.7)log 10 [i?/i? 200 ]+(1.3±0.6) (9) 

where SFR is the average in each radial bin. The fit 
yields xt — 1-2 and is shown by the solid line on the 
upper panel with the la uncertainty given by the grey 
region. The red line sho ws the gradient predict ed by 
the gas starvation model of Book & Benson (2010), nor- 
malized to minimize the \ 2 between the model and the 
observed SFRs. The model predicts the average SFR for 
all non-BCG cluster members at z = 0, so the higher red- 
shift of the cluster galaxies in Figure [7] might explain the 
shift in (SFR) between the observations and the model. 
However, even after we adjust the normalization of the 
model, it remains a poor fit to the data (xt — 3-2), but 
the re-normalized model is marginally consistent with the 
best-fit power law within the large statistical uncertainty 
on the fit. 

Despite the significant uncertainty in the power law 
fit to the SFR-radius relation, Figure [7] clearly demon- 
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strates higher (SFR) toward the outer regions of the 
stacked cluster sample. This mirrors the trend found 
via partial correlation analysis in £16.11 which is signifi- 
cant at >99.9% confidence. A similar trend appears in 
(sSFR) as a function of radius, which yields (sSFR) oc 
(i?/i?2oo) 10±a6 - To interpret Eq. [9] in detail, we need 
a model that better agree s with the observations than 
the lBook fc Bensonl (|201CO model and that accounts for 
projection effects, the distribution of orbits followed by 
cluster members, and the effect of different environmen- 
tal processes. 

The RPS scenario makes at least one clear, qualitative 
prediction that we can use to evaluate its impact without 
a detailed model. Because RPS operates quickly com- 
pared to the cluster crossing time, the radial variation 
in (SFR) should be caused by variations in the fraction 
of SFGs (/sf), and there should be little change in the 
SFRs of individual galaxies. The middle panel of Fig- 
ure [7] shows that fs f declines strongly near the cluster 
center, which is consistent with RPS. We also see lower 
(SFR) among SFGs with R < 0.1 i?2oo compared to 
SFGs with R > 0.1i?2oo (bottom panel). This difference 
is formally significant at > 99.9% confidence, so SFGs 
with R < 0.1i?2oo experience a clear reduction in their 
SFRs as they transition to passive evolution. However, 
a fit to (SFR) sfg versus R/R200 shows no significant 
trend, so the reduction in SFR at R < 0.1i?2oo may be a 
sharp transition rather than a gradual decline. The time 
to cross this region is approximately 100 Myr, which is 
consistent with the hypothesis that the apparent break 
is due to the influence of RPS. The onset of this break, 
however, occurs much closer to the cluster cen ter than is 
usually expected (~ 0.5R™n: iTreu et aT1l2003fl . 

The small region over which RPS has the strongest ef- 
fect might explain the absence of the correlation between 
M* and SFR that is expected under the RPS scenario. 
The most massive galaxies, which are best able to retain 
their gas, are preferentially found at R < 0.1i?2oo- These 
galaxies therefore experience stronger ram pressure than 
typical, less massive cluster members, and their ability 
to retain their gas is cancelled by the increased stripping 
they experience. 

6.4. TIR Luminosity Function 

Another probe of the impact of environment on star 
formation is the TIR luminosity function (LF). The TIR 
LF is sensitive to the frequency of star formation in clus- 
ters and the rapidity with which it is quenched; this 
provides a strong empirical constraint on the types of 
processes that mediate the interaction between individ- 
ual galaxies and the cluster environment. For example, 
IBai et akl (j2009f ) found similar shapes (a and L*) of the 
TIR LFs in the galaxy clusters that they measure com- 
pared to the field galaxy TIR LF. They argue that this 
similarity requires truncation of star formation on short 
timescales compared to the lifetime of star formation in 
individual galaxies. Such rapid transitions are inconsis- 
tent with processes like gas starvation and galaxy harass- 
ment. 

To evaluate the conclusion that RPS dominates the 
evolution of star formation in galaxy cluster members, 
we will examine the TIR LFs of the clusters in our main 
cluster sample. We construct the LF as described in ^5.21 
and the results appear in FigureEJ The main cluster sam- 
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Fig. 8. — Total infrared (TIR) luminosity functions for each of 
the 5 clusters in the main sample. Input galaxies are restricted to 
Mr < — 20 and R < 0.4i?2oo- The Coma cluster LF (blue triangles; 
IBai et al. 2009) is shown for reference. The solid black lines on each 
panel show the field galaxy TIR LF of Perez-Gonzalez et al. ( 2005) 
shifted to the redshift of the cluster and normalized to match the 
observed LF above our nominal completeness limit of 2.1 X 10 10 Lo 
(vertical, dashed line), which is the approximate Ltir expected 
for a normal, spiral galaxy with Mr = 20 mag. This indicates the 
expected distribution of SFRs among field galaxies that enter the 
cluster. The dashed lines mark the completeness limit imposed by 
the requirement Mr < —20. Some variation between individual 
clusters is apparent at the highest luminosities. 

pie contains only 5 clusters, which prevents construction 
of subsamples that have different masses and similar red- 
shifts. Therefore, we cannot reliably identify effects that 
are strongly dependent on cluster mass. 

The dashed vertical line in F igure El marks the ex- 
pected Ltir of a galaxy with the lAssef et al.l (|2010l ) spi- 
ral SED and Mr = —20. This marks the approximate 
TIR completeness limit imposed by the requirement that 
Mr < -20. We call this limit L T l [^ h . This limit is 
representative only, and Figure [8] includes many cluster 
members that have Mr < -20 and Ltir < ^riR h ■ 
This is expected because cluster galaxies have lower 
(sSFR) th a n the field galaxies used to construct the 
lAssef et all (|20100 templates. In fact, 65% of galaxies 
with Mr < — 20 mag and measurable (> 3a) MIR emis- 
sion are less luminous than • This means that 
L t T 1 [ R ish is robust, and the true limit is lower than the 
nominal value established from the spiral galaxy tem- 
plate. To predict the true L l ^^ h , we would need a model 
for the truncation of star formation in clusters, which is 
exactly what we want to measure. To be conservative, 
we restrict our fits to use only bins more luminous than 
^TiR h - Above this limit, we can be confident that the 
weights given by Eq. [5] will correct to the full galaxy 
popul ation. 

Like IBai et all ()2009h . we find that the individual clus- 
ters in Figure [8] have luminosity functions that closely 
resemble the field galaxy LF at their respective redshifts. 
This agreement occurs despite the disagreement between 
the field galaxy LF and the combed cluster sample (Fig- 
ure [9]) . While it is possible that the field galaxy LF at 
the median redshift of the cluster provides a poor es- 
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Fig. 9. — TIR LF of the stac ked cluster samp le (filled squares) 
compared to the LFs derived bv lBai et all (2009!) for Coma (filled 
triangles) and Abell 3266 (open triangles). The solid lines show the 
best fit LFs for the stacked sample (black) and for Coma (blue), 
and the shaded region shows the 68% confidence interval around 
the best fit to the stacked cluster LF. The da shed line shows the 
field galaxy LF of Perez-Gonzalez ct al. (2005), and we have em- 
ployed the prescription of Lc Floc'h ct al. (2005) to correct the 
LF to the median redshift of the combined cluster sample. The 
field galaxy LF is normalized to best match the observed LF above 
the nominal com pleteness limit. The average LF determined by 
IBai et al.l (I2009T) from two high-z clusters (MS 1054-03 and RX 
J0152, red pentagons) is shown for comparison. 

timate of L^ IR among galaxies that fall into clusters, 
any effects of preprocessing in the large-scale structure 
around the cluster should appear in Figure [8] A better 
explanation appears to be that the improved statistics 
in the combined cluster sample illuminate a discrepancy 
that is not visible in the individual clusters due to larger 
observational uncertainties. 

At least some of the variation between the IR LFs ob- 
served in different clusters may be caused by system- 
atic uncertainties in the completeness corrections. We 
can only apply completeness corrections in regions of the 
clusters where we have both spectra and MIR photom- 
etry, so azimuthal asymmetry may be important. The 
LF of the stacked cluster sample averages over several 
selection regions, so it is less subject to this uncertainty. 

We constructed the best-fit Schechter functions to both 
the Coma cluster and the stacked cluster shown in Figure 
OH The Schechter function has the form, 



cf>* 



L 

L ;( 



(10) 



where gives the projected surface density of sources 
at TIR luminosity L, and a and L» are the usual 
Schechter function parameters. We fixed a = —1.41 in 
the fit to the c luster LF, which is the best-fit value for 
the Coma LF (|Bai et al.l I2006T) . ILe Floc'h et~aTl (pOOl 
suggest that the faint end of the LF cannot evolve much 
with redshift for z < 1, so the faint end of the LF in the 
Coma cluster is likely to provide a good estimate of a 
in all galaxy clusters. The best fit to the stacked main 
sample has = (6.6 ± 1.1) x 10 10 L©. 
If clusters rapidly shut off star formation in galax- 



ies that fall in from the field, as IBai et all (|2009f ) con- 
clude, then only galaxies that have recently become clus- 
ter members will have measurable star formation, and 
the TIR LF of a cluster should have L* and a similar 
to the field galaxy LF at the same redshift. There- 
fore, we want to compare L* to the field galaxy LF 
at the media n redshift of the comb ined galaxy sample, 
z med = 0.24. ILe Floc'h et al.l (|2005D found that the field 
galaxy LF evolves as L* oc (1 + z) n , where n = 3.2+g;g. 
iPerez-Gonzalez et al.l (|2005t ) studied the 12/xm LFs of 
field galaxies from z — to z — 3 and found that the 
field galaxy LF at z = 0.1 has L* 12fim = 4.1 ± 1.3 x 1O 9 L 
and a = 1.23 ± 0. 07. We use the prescription of 
iTakeuchi et~aT1 (|2005| ) to convert their L* to a TIR lu- 
minosity, which yields L* TIR = 2.3 x 10 10 Lq at z = 0.1. 
We determine the field galaxy LF at the median redshift 
of the s tacked cluster sample ( z me d = 0.211) with the re- 
sults of ILe Floc'h et all (|2005f ) and fit the normalization 
of the LF to the observed cluster galaxy LF. The result is 
shown as the blue, dashed line in Figure [9l The quality 
of the fit (x 2 = 4.7) is considerably poorer than the fit 
shown by the heavy, black line (x 2 = 0.5), which uses a 
from the Coma cluster and fits for L^ IR and 4>* . While 
the absolute \ 2 values cannot be used to evaluate the 
quality of the fits due to the presence of correlated errors 
in adjacent bins, the Coma-based LF improves the qual- 
ity of the fit by A\ 2 — 4.2 with only 1 additional degree 
o f freedom. 

IBai et ail (|2009l ) found that the luminous ends of the 
TIR LFs of the Coma cluster and A3266 have similar 
shapes and that the L^, IR for these clusters are indistin- 
guishable from the field galaxy LF. The similarity be- 
tween the stacked cluster LF and the redshifted field 
galaxy LF for L tje > 4 x 10 10 L W is consistent with 
the conclusions of IBai et al.l ()2009| ). They argue that the 
similar LF shapes in clusters and in the field suggests 
that gas starvation is not a plausible mechanism to end 
star formation among cluster member galaxies. Because 
gas starvation operates slowly (~ Gyr timescales), they 
conclude, it should produce many galaxies in the tran- 
sition phase between SFGs and passive evolution. They 
found no such transition population. However, the com- 
bined cluster LF shown in Figure [9] displays a 4tr deficit of 
galaxies with moderate SFRs (SFR « 5 M yr -1 ) com- 
pared to the expectation from the field LF. This largely 
drives our conclusion that an unmodified field galaxy LF 
does not provide a good description of the SFG popu- 
lation in the merged cluster sample. The disagreement 
between the field and cluster galaxy LFs could indicate 
the presence of a transition population. It is possible 
that the discrepancy results instead from some selection 
effect not accounted for in our completeness estimates. 
The most obvious culprit for such an effect is some resid- 
ual dependence of spectroscopic member identification 
on color that we have not been able to identify from the 
data. 

As a test for radial gradients in the population of tran- 
sition SFGs, we binned the galaxies in the main clus- 
ter sample into three radial bins with equal numbers of 
galaxies. The TIR LFs for the radial subsamples are 
shown in Figure 1101 The L^ IR increases slightly from 
the innermost to outermost radial bins, but this increase 
is not statistically significant. This marginal decrease in 



12 



Atlee et al. 




10 10 10" 

L(8/im-1000Mm) [L s ] 



Fig. 10. — TIR LF divided into radial bins with equal numbers 
of galaxies. Each panel is labeled with the range of radii that 
contribute to the LF, where r = R/Raoo- The solid lines show 
the best-fit Schechter function to the LF in each radial bin, and 
the shaded regions show the uncertainties on the fit. The fits are 
constructed from the bin more luminous than the completeness 
limit (dash e d lin e), which is set by the expected Ltir for the 
lAssef et alj (2010) spiral galaxy template at Mr = —20 mag. The 
red, dashed lines on each panel show the best-fit Lyj„, and its la 
confidence interval is given by the shaded region. The three LFs 
hint at an increase in L%,j„ analogous to the increase in (SFR) sfg 
seen in Figure [7] but this increase is not statistically significant. 

L^ir m the in nermost radial bi n is qualitatively similar 
to the results of lBai et al.l (|2009f) , who also examined the 
radial dependence of the TIR LF with a very similar bin- 
ning scheme. However, the fractional change in Lrpjj^ in 
the Coma cluster is much larger than observed in Figure 
[T0l The increase in (SFR) sfg seen in Figure [7] in the 
same radial bins does not appear to affect L^ IR . This 
indicates that the change in (SFR) sfg is driven by the 
frequency of galaxies with SFR s» 3 Mq yr' 1 . 

6.5. Substructure and Preprocessing 

Von der Linden et al. (2010) found a trend toward in- 
creased (SFR) at larger R/R200 that extended out to 
at least 2i? 2 oo- They concluded that preprocessing in 
groups contributes significantly to the SFR-density rela- 
tion. Our observations do not extend past R = 0.4i? 2 oo, 
so it is impossible to measure preprocessing directly, but 
A3128 shows significant substructure, so we can compare 
it to the smooth clusters in the sample to probe how 
the presence of substructure influences SFRs in clusters. 
This allows us to indirectly test the impact of group scale 
environments on SFGs, because coherent substructures 
in clusters should correspond to recently-accreted groups. 

Before we can compare the integrated SFRs in dif- 
ferent clusters, we must first correct for the different 
numbers of galaxies in each cluster. The sSFR natu- 
rally accounts for this variation, and it is therefore a 
better parameter to compare integrated SFRs between 
clusters. We employ a method analogous to Eq. [5] to 
calculate the (sSFR) and compare A3128, which shows 
significant substructure, to the other clusters in the main 
sample. We find (sSFR) = 9.0t^ x 10~ 12 yr" 1 and 
(sSFR) = 3-l±g;l x 10- 11 yr- 1 in A3128 and in clus- 



ters without substructure, respectively. If we correct 
(sSFR) of A3128 to the mean redshift of the other clus- 
ters (z = 0.241), we find (sSFR) = 1.5t° A 4 x lO" 11 yr' 1 . 
This is still lower than the average of the clusters without 
substructure at > 99.9% significance. 

The difference between A3128 and the other clusters 
might be a result of the structure in A3128, or A3128 
might simply have an unusually low (sSFR) for its red- 
shift. In the latter case, the observed difference would 
be a result of cosmic variance. We compared A3128 with 
the 4 individual clusters without substructure, and we 
found that A3128 has higher redshift-corrected (sSFR) 
than MS1008, which is approximately 50% more mas- 
sive than A3128. However, the typical dispersion in 
Jsf amon g nearby cluste r s with a > 800 km s' 1 is 
~ 0.1 dex iPoggianti et "all (|2006() . so to explain the ob- 
served deficit of (sSFR) in A3128 as cosmic variance 
would require a ~ 3<r excursion. While cosmic variance 
provides a marginally plausible explanation for the deficit 
of (sSFR) in A3128, the presence of substructure ap- 
pears to be the more likely cause. If the observed differ- 
ence arose from groups that have recently fallen into the 
cluster, the excess sSFR in clusters without substructure 
would imply that the "average" group member is likely to 
have experienced preprocessing. This result may be ab- 
sent from the partial correlation results (Tabled]) because 
only ~ 10% of cluster members have ever been part of a 
large group (|Berrier et a l. 2009), and only former group 
members that have been accreted more recently than the 
e-folding timescale for star formation will show evidence 
of preprocessing. 

The conclusions drawn from the measured SFRs of 
cluster members and from the integrated cluster prop- 
erties depend strongly on the methods used to identify 
and correct for AGN. In Paper I, we noted that the IR 
and X-ray AGN selection techniques identify quite dif- 
ferent samples. If we relied only on X-ray based AGN 
selection, as some authors do, the MIR luminosity con- 
tributed by unidentified AGN would lead to an overes- 
timate in the integrated SFR of the cluster. For exam- 
ple, in A1689 we would overestimate the total SFR by 
20%. Applied to all clusters simultaneously, this alter- 
native method of AGN correction results in an inferred 
(sSFR) = 7.2tj 3 x 10~ n yr- 1 among the clusters with- 
out measurable substructure but no measurable change 
in A3128. In this example, uncorrected AGN contamina- 
tion would dominate the observed difference in (sSFR), 
and we would over-estimate the impact of preprocessing 
in the group environment. 

6.6. MIR Butcher- Oemler Effect 

The relative importance of gas starvation and RPS 
is also probed by the evolution in (SFR) as a func- 
tion of cosmic time. The classic example of this i s 
the Butcher - Oemle r effect ()Butcher fc Oemler! 11978ft . 
lHaines et al.l (|2009l ) constructed an analogous measure- 
ment with SFRs measured via vL v {1\[im) among the 
LoCuSS cluster galaxies. They employed a SFR thresh- 
old of 8.6 Mq y^ 1 , and they found that fsF oc (1 + z) n 
with n = 5.7j^"g. Figure [TT] shows their fit to fsF among 
the LoCuSS clusters as a function of redshift. The fsF 
values for our cluster s and for a higher red shift cluster 
sample measured by iSaintonge et al.l (|2008l ) are super- 
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Fig. 11.— Integrated fractions of SFGs (f SF , SFR > 
8.6M0|/r _1 ) in all 8 clusters as a function of redshift. Filled red 
squares mark the clusters in our sample for which we successfully 
measure /sf, an d red arrows mark the clusters for which we can 
produce only upper limits. Open red squares mark the f$F that 
would be inferred from X-ray only AGN identifications. The filled 
and open squares overlap for all clusters except AC11 4. Open grey 
triang les indicate the LoCuSS clusters as reported bv lHaines et al.l 
1200 91), and open blac k pentagons mark the clusters measured by 



ISaintonge et all H2008IV The soli d line indicates the best-fit f$F - 
relation from Haines et al. (2009), and the shaded region shows the 



la confidence interval for their fit (fsF c< (1 + z) 
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imposed. The 8 clusters in our sample, shown as the 
red triangles in Figure HT| are clearly consistent with 
the IHaines et a l. (2009) result within the uncertainties. 
However, the fit to the LoCuSS cluster s syste matically 
overpredicts Jsf in the ISaintonge et all (|2008l ) clusters, 
despite the lower SFR threshold (5 M Q yr~ 1 ) used by 
ISaintonge et~aTl (|2008l ). 

In £16.51 we considered the impact of X-ray only 
AGN identification on the inferred (sSFR). This be- 
comes a more important consideration at high-z, because 
the frequency of lum inous AGNs increases dramatically 
([Martini et al.l [2(501 . Fi gure Qj] includes two points for 
each cluster in our sample. One shows Jsf with the IR 
AGN selection included (filled triangles), and the other 
shows fs f that we would measure if we only knew about 
the X-ray selected AGNs (open triangles). The Jsf in- 
ferred from the X-ray only selection in AC 114 differs by 
1.6<7 from the result when the full AGN sample is consid- 
ered. This illustrates the contamination that X-ray only 
AGN identification can introduce to integrated SFRs. 
This contamination becomes more severe, and appears 
in other clusters, for SF R thresholds l ess th an the fairly 
high value employed by IHaines et al.l ()2009l ) . 

7. DISCUSSION 

In £16.11 and £|6.3l we examined correlations between en- 
vironment, SFR and M*. We found a strong correlation 
of SFR with i?/i?200, and we found evidence for a tran- 
sition population of low-SFR galaxies near the cluster 
center. We interpret this population as evidence that 
galaxies in this region experience a rapid reduction in 
SFRs that initiates their transition from SFGs to passive 
galaxies. This interpretation was supported by a possi- 



ble trend toward larger L^, IR farther out in the cluster 
f £16.4[) . We also found evidence for a concentration of 
massive galaxies near the cluster center f £16.2[) . 

In this section, we consider the results of i i6.3ll6.5l in 
more detail and interpret them in the context of two com- 
peting mechanisms to end star formation in cluster galax- 
ies: RPS and gas starvation f £|7.1[) . We also briefly dis- 
cuss the additional information that the Butcher-Oemler 
Effect can provide abou t the impact of the cluster envi- 
ronment on SFGs (JL2]). 

7.1. Star Formation in Clusters 

In £J51 we examined several diagnostics for the impact 
of the cluster environment on star formation. These in- 
clude partial correlation analysis, (SFR) and (sSFR) 
versus radius, and an examination of TIR LFs. The 
partial correlation results informed much of our subse- 
quent analysis. One important result was the absence of 
a correlation between SFR and £in once we control for 
R/R2Q0- This implies that interactions between individ- 
ual galaxies have only a limited impact on star formation 
in cluster members. We conclude that the SFRs of clus- 
ter members are controlled by hydrodynamic interactions 
between galaxies and the ICM. 

There is disagreement in the literature concerning the 
importance of different mechanism s to shut down star 
formation in clusters. For example, iSimard et al.l (2009) 
determined that evolution in cluster SFRs is controlled 
by galaxy-galaxy interactions because the growth in the 
fractions of early-type and passive galaxies track one 
another very closely in their sample. This contrasts 
sharply with our results in £16.11 which suggest that in- 
teractions with the ICM are the dominant factor. Of 
the hydrodynamic processes commonly considered (e.g. 
RPS and gas starvation), only RP S has been directly 
observed to work in nea rby clusters (|Kennev et al.ll2004l : 
iSivanandam et al1l2010l ). We therefore ask whether our 
observations are consistent with RPS alone or if an addi- 
ti onal mechanism is required to explain the observations. 

iTreu et all (|2003| ) determined that RPS works effec- 
tively for Milky Way-like galaxies in a cluster with 
M vir = 8 x 10 14 M Q when R < 0.5i? 2 oo- Our sample 
is restricted to projected R < 0.4i? 2 oo, and their clus- 
ter mass is similar to the typical cluster in our sample 
(M clust w 5 x 10 14 M Q ), so RPS should act efficiently 
on most galaxies in our sample. Nevertheless, some of 
the obvious signatures of RPS do not appear: We found 
no residual correlation of SFR with M* at fixed -R/-R200, 
which is contrary to the prediction that RPS should af- 
fect low-mass galaxies more strongly. Furthermore, we 
found in £16.41 that the TIR LF of cluster galaxies is sig- 
nificantly better fit by a Schechter function with variable 
^tir an d with the Coma cluster's faint-end slope than 
by a redshift-corrected field galaxy LF. This is inconsis- 
tent with the simplest prediction of RPS, which would 
suggest that the SFG population in cl usters should be 
well-d escribed by the field galaxy LF (|Bai et al.l 120061 
2009). To interpret this disagreement, we need to know 
whether it results from different faint-end slopes between 
the cluster population and the field, different L^ IR , or 
some combination of the two. This level of detail is 
impossible given the completeness limit in the present 
sample, so the implications of the poor agreement be- 
tween the cluster and field SFG populations are ambigu- 
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ous. Finally, the best-fit power law to the SFR-radius 
relation is consistent — within very lar ge observational 
uncert ainties — with the predictions of Boo k fc Benson! 

However, we also found evidence that RPS contributes 
significantly to the decline of SFRs among cluster galax- 
ies. Figure [7J shows that the dependence of SFR on 
R/R200 is driven largely by a decline in Jsf toward the 
cluster center, and only a small fraction of the depen- 
dence is driven by decline in the SFRs of individual SFGs. 
This is more consistent with RPS than with gas starva- 
tion. Furthermore, the residual decline in (SFR) sfg 
with R/R200 does not occur smoothly, but appears to 
set in rapidly at R 0.1i?2oo- The crossing time for this 
region (~ 100 Myr) is consistent with RPS as the origin 
of the transition, but the small size of this region com- 
pared to the radius where RPS is theoretically expected 
to be important (0.5 R200) is surprising. This explana- 
tion is also difficult to reconcile with the decline in fsF, 
which we attribute to RPS, across all radii. 

To interpret Figures [Jj and 1 1 01 we must consider pro- 
jection effects and the influence of galaxy "backsplash" . 
Projection effects will cause some galaxies at intrinsically 
large R/R200 to appear at small radii when the cluster 
is projected onto the plane of the sky. If we assume an 
R~ 2 density p rofile for galaxies and fsF v s. R/R200 as 
determined bv lvon der Linden et al.l (| 20101 ). only ~ 30% 
of SFGs with projected R < 0.3i? 2 oo actually fall within 
that region. This re p resen ts an upper limit, since the 
Ivon der Linden et al.l ([2010) result also includes projec- 
tion effects. Therefore, we cannot assume that SFGs at 
small R/R200 actually reside in the high density regions 
where RPS is most important. This suggests that SFGs 
that physically reside inside O.I-R200 have their SFRs re- 
duced more drastically than implied by Figure [71 

"Backsplash" refers to galaxies on nearly radial orbits 
that pass through the dense central region of the cluster 
and return to large i?/i?2oo- This effect can make radial 
gradients like the ones shown in Figure [Jj particularly 
difficult to interpret, because even galaxies presently at 
large radii may have passe d near the cluster center in 
the past. IGill et aT] (|2005f ) use N-body simulations to 
find that 50% of galaxies between l-2i?2oo are backsplash 
galaxies, and 90% of thes e have be e n insi de 0.5i?2oo at 
some point in the past. Pimbblet (201l|) use mixture 
modeling with observations of real clusters from SDSS to 
infer that 60 ± 6% of galaxies at i?/i?200 = 0.3 are part 
of the backsplash population, so they were even deeper 
into the dense central region of the cluster at some point 
in the past. Galaxies on radial orbits in a cluster with 
an R~ 2 mass profile and a core radius of 0.05i?2oo spend 
only 4% of their time inside O.I-R200 and only 32% of their 
time inside 0.5i?2oo- This model also implies that approx- 
imately half of the galaxies in a cluster with an age of 8 
Gyr have previously passed through the cluster center. 
As a result, the dynamic nature of the cluster populations 
smears the dependence of star formation on projected ra- 
dius relative to the underlying, three-dimensional trends. 
This suggests that the projected trends shown in Figure 
UJare lower limits to the true, three-dimensional trends. 

The effects of projection are potentially just as impor- 
tant. We adopt the simple mass profile described above, 
which implies that ~ 40% of galaxies with projected 
R < O.I-R200 actually reside outside 0.1i? 2 oo- However, 



the trend in fs f versus radius (Figure [JJ implies that at 
least 67% of SFGs that appear with R < 0.1i?2oo in pro- 
jection are likely to have three-dimensional radii outside 
this limit, and 20% of SFGs with projected R < 0.4i? 2 oo 
will lie outside this radius. The steep dependence of fsF 
on radius therefore works to mitigate any radial trends 
among the SFG population, which would lead to the 
small dependence of (SFR)sfg on R/R 2 oo- This sug- 
gests that the sharp drop in (SFR) sfg at projected ra- 
dius O.I-R200 would be stronger if measured relative to 
physical radius. 

An ideal way to account for both projection and back- 
splash in our observables is to compare our results to 
models that include these effects. This approach allows 
more reliable conclusio ns than simple, ad hoc arguments. 
iBook fe Benson! ((2010) developed a model for the re- 
moval of hot gas from galaxies by the ICM, which is the 
physical mechanism that drives gas starvation. Their 
"shocks" model predicts that galaxies that experience 
this process should show SFR oc (i?/i?20o)~ ' 6 between 
0.1-0.4i?2oo (their Figure 3). This is consistent with 
our results in gEJ ((SFR) cx (i?/i?20o) l ;3±0 - 7 )- How- 
ever, the model yields overall poor agreement with the 
data, even after we fit the normalization of the model 
to best m atch the observa t ions. This contrasts with the 
results of IBook fc Bensonl (|2010f ). who found that their 
model agrees well with the SF R-radius relation m easured 
among the CNOC clusters (jBalogh et al.l I2000D . The 
resolution of this conflict will require additional obser- 
vations and more sophisticated theoretical models that 
provide predictions for competing processes. 

The large statistical uncertainties on the measured 
SFR versus R/R200 preclude detailed comparisons be- 
tween our observations and a model for either RPS or 
gas starvation, so we must rely on qualitative arguments. 
Figure [9] demonstrates that the field galaxy TIR LF at 
the median redshift of the combined cluster galaxy sam- 
ple provides a poor match to the observed TIR LF, and 
this discrepancy is most pronounced among galaxies with 
the lowest SFRs. The bottom panel of Figure [7J also 
shows a sudden decrease in the SFRs of SFGs with pro- 
jected radii R < 0.1i?2oo- The crossing time for this 
region is ~ 200 Myr, which suggests that RPS is respon- 
sible for this reduction. However, only 5% of SFGs in 
the cluster have R < O.I-R200, so this reduction in the 
SFRs of SFGs near the cluster center cannot account for 
the discrepancy between the TIR LFs of cluster and field 
galaxies. There must also be an effect at larger radii that 
is not apparent in the present sample. The entire sample 
has R < 0.5i?2oo, so either RPS or gas starvation could 
plausibly cause the disagree ment. 

In contrast to our results, iBai et all (|2009f) found that 
the TIR LFs of many clusters are consistent with one 
another and with the field LF. They inferred that cluster 
galaxies only rarely occupy a transition phase between 
SFRs characteristic of field galaxies and complete pas- 
sivity. From this, they determined that star formation 
in cluster galaxies must be truncated on short timescales 
compared to th e lifetime of the cluster. Both we and 
IBai et ail (|2006f ) find ~ la variations in the shape of 
the TIR LF in different R/R200 bins. The observed de- 
creases in L^p IRl while not statistically significant, are 
consistent with the decline in (SFR) sfg for R < O.I-R200 
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(Figure [7]). The latter result indicates that the high 
ICM density near cluster centers reduces SFRs in indi- 
vidual SFGs. The smooth decline in fgp as a function 
of R/R200 suggests that this process eventually results 
in the end of star formation in these galaxies. Projection 
effects and backsplash will both influence the observed 
trends. Projection would cause us to mistake galaxies at 
large R/R200 for galaxies near the cluster center, while 
backsplash would move galaxies that had been processed 
near the cluster center back to the outskirts of the cluster. 
Both effects would cause the projected trends to appear 
weaker than the true, three-dimensional variations in the 
cluster. This suggests that the observed radial variations 
are real, and the trends with projected radius likely un- 
derestimate the intrinsic, three-dimensional trends. 

We find lower (SFR)sfg inside O.I-R200 compared to 
outside this radius. We also see hints of this change 
in the lower L^ IR near the cluster center compared to 
further out. The variation in the properties of SFGs 
implied by these measurements indicates a substantial 
change in the SFGs very close to the cluster center. The 
crossing time for the sphere with radius O.I-R200 is less 
than 200 Myr, which strongly favors RPS as an explana- 
tion. However, we also see indications for a deficit in the 
number of low-SFR cluster galaxies relative to the field 
population. Because only ~ 5% of SFGs have projected 
R < 0.1i?200i galaxies outside 0.1i?2oo must dominate the 
under-abundance of galaxies with SFR ss 3Mq yr^ 1 . 

We can determine whether RPS or gas starvation are 
more likely to be responsible for the deficit from the time 
required to make the transition. We find that 66% of 
SFGs with Mr < -20 mag have L TIR < Lf™ sh , where 
^ t TiR h 1S the luminosity expected for a typical field spi- 
ral with M R = -20. If 50% of field SFGs with the same 
M R distribution had L T ir < Lffff* 1 , then 16% of clus- 
ter SFGs would be in transition. Combined with the 
gas c onsumption timesc ale of a typical spiral galaxy (2.4 
Gyr; iBigiel et all 120111) , this implies a transition time 
of ~ 400 Myr. This timescale is approximately twice 
the dynamical time of an ordinary spiral galaxy. This 
timescale favors RPS as the dominant mechanism among 
the full cluster sample, because RPS implies that galax- 
ies should remain in a transition phase for approximately 
their dynamical time while the cold ISM is stripped. 
However, the assumption that 50% of field SFGs in an 
M^-matched sample would have Ltir < ^tir 1S arbi- 
trary. A comparison of the SFR-Affj relations in clusters 
and in the field is required to measure the transition time 
more precisely. 

If star formation in most cluster galaxies ends as a 
result of RPS, post-starburst galaxies should be more 
frequent in clusters than in the field. This is a robust 
prediction of any scenario that results in a rapid transi- 
tion of SFGs to passive evolution. Galaxies with K+A 
spectra, which are usually associated with post-starburst 
populations, should remain visible for ~ 100 Myr to 
1 Gyr. This is short compared to the cluster crossing 
time, so a large population of K+A galaxies relative to 
SFGs would be str ong evidence that RPS plays an impor- 
tant role. Instead, lYan et afl (|2009f ) report that galaxies 
with K+A spectra are less common in overdense envi- 
ronments like clusters than in the field at z rs 0.1. They 
suggest that K+A galaxies appear at constant absolute 



density, and th at this density corres ponds to the group 
scale at z 0. Dressie r" et al.l (|1999f ) instead found that 
the fraction of K+A galaxies is much higher in c lusters 
than in the field, and lvon der Linden et al.l (|2010f ) found 
no dependence of the ratio of Nk+a/Nsf on i?/i?200- 
The different metho ds used by the s e aut hors to select 
their K+A samples — [Dressier et al.l (119991) rely on [On] 
to exclude SFGs, wh i le lYan et al.l (120091 ) use H/3, and 
Ivon der Linden et al. I (l2010ir select galaxies with excess 
Balmer line absorption from their principal component 
analysis — may account for the apparent contradictions 
in these observational results. 

Groups that have recently fallen into a cluster might 
appear as an excess in t he substructure parameter 
(|Dressler fc Shectmanl lT988). In our sample only A3125 
has a mass comparable to galaxy groups, and we have 
not considered that cluster in our analysis. Therefore, 
we cannot directly constrain the mechanism that drives 
SFR evolution in group members. However, we find that 
(sSFR) is higher among clusters with no substructure 
than in A3128, which is the only member of our main 
sample with significant substructure. This could indi- 
cate that galaxies that have recently been part of groups 
have experienced preprocessing, but it could also arise 
from cos mic variance. In a rec e nt stu dy of SDSS galaxy 
clusters, Ivon der Linden et al.l (|2010l) found a trend of 
SFR with radius that extended to 2i?200- They con- 
cluded that preprocessing of galaxies before they become 
cluster members is likely to contribute significantly to 
the SFR-radius relation. Additional observations are re- 
quired to determine if galaxy groups are responsible for 
the preprocessing or if other processes are required. 

7.2. Evolution 

In Ej6.6l we suggested that the evolution of star forma- 
tion in clusters is sensitive to the mechanism(s) responsi- 
ble for the appearance of the z = SFR-density relation. 
In particular, the rate of evolution of Jsf is sensitive 
to the operation of the cluster environment on recently 
accreted field galaxies. Figure [TT] shows the fraction of 
SFGs in clusters as a function of redshift since z sa 0.8 for 
th e clusters in our s ample (red ) compared to the sa mples 
of iSaintonge et al.l (|2008l ) and (Haines et al.l (|2009l ). 

We use measurements of Jsf versus z to estimate the 
time required for Jsf in clusters to decline by a fac- 
tor of e compared t o coeva l field galaxies. The best- fit 
to the IHaines et all (|2009D galaxies (f SF oc (1 + z) m , 

m = 5.7^'g) is shown as the black line in Figure 
TT1 and it agre e s well with the clusters in our sample. 
Le Floc'h et all (|2005l ) report that the field galaxy LF 
evolves as L^ IR (z) oc (1 + z) n , where n = 3.2^0.2- Ap- 
proximately 70% of cluster member galaxies at z = 
never had a massive companion before t hey entered the 
cluster environment (|Berrier "eTaIll2009l) . so we can as- 
sume that galaxies that fall into the cluster have the same 
LF as field galaxies. The threshold we use to identify 
SFGs (SFR > 8.6 M Q yr" 1 ) is larger than the SFR that 
corresponds to L^ IR (4 M Q yr" 1 ), so we assume that 
fsF among field galaxies has the same redshift depen- 
dence as L^jr- With this assumption, we can measure 
the relative change in Jsf as a function of redshift and 
determine how the cluster environment induces SFGs to 
turn passive. The ratio of fsF,dust to fsFjieid has un- 
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dergone approximately 1.7 ± 1.2 e- foldings since 2 = 1. 
The elapsed time over this redshift interval is 7.7 Gyr, so 
the e-folding time for fsF,clust/ fsFjieid is 4.6± 1 ° 8 6 Gyr. 

The e-folding time of fsF does not correspond directly 
to the truncation time for star formation in individual 
cluster members. New SFGs constantly fall into the clus- 
ter from the field, and this results in a longer timescale for 
f SF,dust/ fsf, field to decline than for SFRs to decline in 
individual galaxies. The rate at which SFGs fall into the 
cluster combines with the timescale for the conversion of 
individual SFGs to passive evolution to determine how 
rapidly fsF,dust/fsF,field changes. The timescale for 
this change is long compared to t he gas exhaustion time 
in typical spiral galaxies (2.4 Gvr. lBigiel et al.|[201 lh . but 
this does not necessarily indicate that the timescale for 
the evolution of individual SFGs is similarly long. If the 
timescale for evolution of individual SFGs is indeed long, 
it would favor gas starvation over RPS as the primary 
mechanism to end star formation in cluster galaxies. 

In addition to the degeneracy between changes in 
rates of infall and the timescale for individual SFGs 
to stop forming stars, the measured timescale for 
fsF,dustl fsFjieid to evolve includes sign i ficant obser- 
vational uncertainty. The lHaines et al.l (|2009f ) best- 
fit overpredicts fsF among the high- z clu sters, despite 
the h igher SFR threshold employed by lHaines et al.l 

(2009) , so it underestimates the timescale over which 
fsF,dustl fsF, field evolves. Additional observations are 
required to correct this bias. A measurement of fsF 
versus z with a longer redshift baseline and consistent 
identification of SFGs will appear in our next paper. 

8. SUMMARY AND CONCLUSIONS 

We have used visible to MIR observations of 8 low- 
z galaxy clusters to constrain the impact of the cluster 
environment on star formation. We examined the rela- 
tionship between star formation and environment among 
cluster members and found a strong correlation, with 
(SFR) oc (i?/i?20o) 1 ' 3±0 ' 7 J and this si mple power law 
provid es a good match to the data. IBook fc Benson! 

(2010) model the impact of gas starvation on star for- 
mation. Their model is marginally consistent with the 
power-law that best fits the observed SFR-radius rela- 
tion, but it is a poor fit to the data themselves. The 
(SFR)-R/R2oq relation is dominated by a decline in 
the fraction of SFGs toward the cluster center, but 
we also find lower (SFR) among SFGs with projected 
R < 0.1i?2oo- The dominance of the decline in fgp and 
the short crossing time of a sphere with radius O.I-R200 
both suggest that RPS contributes significantly to the 
observed trend in SFR with R/R 2 oo- The TIR LFs hint 
at a shift toward lower Lfp IR near the cluster center, 
which is consistent with the observed decline in (SFR), 
but this shift is not statistically significant. Projections 
effects and backsplash both work to weaken the observed 
trends relative to the intrinsic variation in three dimen- 
sions, which can hide the steep gradients that would be 
expected from RPS. 

We also examined the relationship between i?/i?2oo 
and stellar mass in cluster galaxies. We found that 
galaxies with R < O.I-R200 show larger (M*) than galax- 
ies farther out in the cluster, even after we have elim- 
inated BCGs from our sample. This excess is signifi- 
cant at ~ 3.5cr, and projection effects are also expected 



to weaken the observed trend relative to the intrinsic, 
three-dimensional variation, so we conclude that it is ro- 
bust. Von der Linden et al. (2010) found no such excess 
once they had removed BCGs, so our result conflicts with 
theirs. This difference may result from the SDSS fiber 
collisions. Our sample is limited to galaxies more lumi- 
nous than the SPS S r-band magnitude limit at the me- 
dian redshift of the Ivon der Linden et al.l (|2010f ) cluster 
sample, so our sample is on average more massive than 
theirs. However, the timescale for dynamical friction to 
affect cluster members is much longer than the Hubble 
time. This suggests that cluster galaxies undergo mass 
segregation via virial relaxation, analogous to the mass 
segregation exhibited by some Galactic globular clusters. 

We measured the fractions of SFGs in our cluster sam- 
ple as a function of redshift, and we found that these 
fractions are consis tent with the measurements made by 
lHaines et al.l (|2009l ) for the LoCuSS clusters. However, 
incomplete AGN subtraction can introduce significant 
contamination to the integrated star formation in galaxy 
clusters. For example, we found that eliminating only 
X-ray AGNs from the sample prior to calculation of fs f 
results in a ~ la excess in AC 114 in Figure ITTl The con- 
sequences are both more significant and more widespread 
for lower SFR thresholds. This can bias measurements 
of star formation as a function of redshift, since the AGN 
contribu tion is expected to b e more significant at higher 
redshift (Mar tini et al.l [2009h ■ With a long enough red- 
shift baseline, evolution in fgp with cosmic time can 
probe the timescale for the end of star formation in clus- 
ter galaxies. Present observations favor gas starvation 
over RPS, but these include important systematic un- 
certainties. 

On balance, our measurements are most consistent 
with RPS as the primary mechanism to reduce star for- 
mation in cluster galaxies. The SFR-radius relation 
agrees better with the predictions of RPS than with gas 
starvation over the range o f radii that we study. This 
supports the conclusions of iBai et al.l (|2009f ) . who stud- 
ied star formation in clusters over a similar range of 
radii and determined that RPS dominates the reduction 
of SFR among cluster galaxies. The observation that 
SFR does not correlate strongly with M* among SFGs 
in clusters can be explained by the concentration of mas- 
sive galaxies near the cluster cent er, where RPS oper- 
ates most efficiently. Furthermore, IBai et ail (|2009T ) also 
found that the luminous end of the TIR LF does not 
vary significantly between z = clusters, and we sim- 
ilarly find that the redshift- appropriate field galaxy LF 
provides a good match to the observed TIR LF in each 
cluster. While this is also consistent with RPS as the 
primary mechanism to end star formation among cluster 
galaxies, the disagreement between the field galaxy LF 
and the stacked cluster sample suggests that there are 
small deviations between the individual clusters and the 
field galaxy LF that we overlook due to limited precision 
in the measured LFs. Such deviations would be more 
consistent with gas starvation than with RPS. Addition- 
ally, the long timescale we infer for the evolution of fs f 
conflicts with the conclusion that RPS ends star forma- 
tion among cluster members and is more consistent with 
gas star vation. This evidenc e ag rees with the conclu- 
sions of IVerdugo et al.l ()2008l ) and Ivon der Linden et al.l 
( 2010), who found independent evidence in favor of gas 
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starvation. A measurement of the rate of change in fgp 
as a function of redshift can provide an additional line of 
evidence to help resolve this disagreement. Present re- 
sults favor a long timescale, but these include significant 
systematic uncertainties. A measurement with a single 
sample of uniformly analyzed clusters will be the subject 
of our next paper. 
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APPENDIX 
ASYMMETRIC DISTRIBUTION FUNCTION 

In order to apply the method described in Jj5]to construct a luminosity function, we must smoothly distribute the 
weight of a galaxy across the specified luminosity bins. The method described in £15.11 to calculate Ltir produces 
asymmetric uncertainties, Ltir = H-a]' > so we nee d an asymmetric probability density function (PDF) to distribute 
weights correctly. This PDF must reduce to the Normal distribution in the case when the upper and lower luminosity 
uncertainties are equal (i.e. Gaussian errors). Here, we describe a piecewise smooth function that satisfies these 
requirements. 

First, wc define an effective dispersion a e = ^Jo\a u , where a u and 07 are the upper and lower uncertainties on Ltir, 
respectively. We then define an alternative dispersion, cr(L), which describes the instantaneous shape of the PDF at 
a luminosity L, 



a(L) 



01 



&e + {o~ u - cr c 



\L-jA 



Al-iA 



IF L < fi - ai 

IF fi — ai < L < fi 

IF L = fx 

IF fj, < L < fx + a u 
IF fj, + o- u <L 



(AI) 



where fi is the best estimate Ltir', o~ u and 07 are the upper and lower uncertainties on fi, respectively. cr(L) smoothly 
connects the low-L and high-L tails of the desired distribution function. Given a(L), we can calculate the probability 
density for a galaxy with measured luminosity fi at L. This probability density is given by, 



f(L,fi,a u ,ai) 



1 



-(p-L) 2 /2a 2 (L) 



2ir cr(L) 



(A2) 



where c(L) is given by Eq. I AH 

The PDF described by Eqs. IA1I and IA2I approaches Gaussian at the high- and low-L extremes, with dispersions a u 
and o~i, respectively. It also smoothly connects these two limiting cases, integrates to unity, and has dispersion equal 
to the geometric mean of o~i and o~ u at the nominal luminosity. It therefore gives a PDF for the luminosity of a given 
galaxy that satisfies our requirements and that is consistent with the available information about Ltir- 
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TABLE 1 

Spectroscopic Completeness 





m R (Vega) 


R/R200 


fci , spec 


°7 


a3128 


14.33 


0.11 


0.669 


0.016 




15.76 


0.10 


0.682 


0.042 




16.32 


0.07 


1.020 


0.065 




16.89 


0.06 


0.532 


0.057 


a3125 


14.60 


0.17 


0.000 


0.000 




15.92 


0.22 


0.000 


0.000 




16.29 


0.05 


1.000 


0.000 




16.77 


0.16 


0.000 


0.000 


a2104 


17.28 


0.05 


0.719 


0.015 




18.32 


0.07 


0.249 


0.007 




18.88 


0.11 


0.222 


0.008 




19.31 


0.11 


0.068 


0.003 


al689 


17.02 


0.03 


0.626 


0.014 




18.46 


0.08 


0.409 


0.014 




19.07 


0.09 


0.217 


0.015 




19.65 


0.12 


0.411 


0.033 


msl008 


19.10 


0.08 


0.693 


0.054 




19.79 


0.04 


0.293 


0.013 




20.35 


0.08 


0.225 


0.014 




20.91 


0.11 


0.084 


0.010 


ac!14 


19.00 


0.10 


0.744 


0.069 




19.84 


0.06 


0.925 


0.072 




20.37 


0.13 


0.399 


0.050 




20.83 


0.14 


0.913 


0.353 



Note. — A sample of spectroscopic completeness 
measurements as described in £|4.1| The complete tabic 
is available from the electronic edition of the journal. A 
brief sample is shown here for guidance regarding form 
and content. 



Atlee et al. 



TABLE 2 

MIR Completeness 





R/R200 


f v (8/j.m) [Jy] 




M24AH71) [Jy] 


C24 fim 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



A3128 



A3125 



A2104 



A1689 



MS1008 



AC114 






.06 


2. 


.65 


x 


10" 


-4 





00+" 


15 
.00 


1 


.70 


x 


10" 


-4 


0. 


oo+"-^ 

— 0.00 





.06 


8. 


.44 


X 


10" 


-4 





.51+" 

— 


18 
17 


5 


.41 


X 


10" 


-4 





■oo+!!'^ 

— 0.00 





.06 


2 


.69 


X 


10" 


-3 





.55+!! 

— 


14 

20 


1 


.61 


X 


10" 


-3 





01 +0.16 
— 0.17 





.06 


8. 


.57 


X 


10" 


-3 





55+"' 
J — 


14 
20 


5 


.44 


X 


10" 


-3 


1 


00+"?" 

— 0.19 





.06 


2 


.66 


x 


10" 


-2 





55+" 


14 
20 


1 


. 74 


x 


10" 


-2 


1 


00+"?" 

W — 0.19 





.06 


8. 


.48 


X 


10" 


-2 





55+"' 

J — 


14 
20 


5 


.51 


X 


10" 


-2 


1 


00+"?" 

— 0.19 





.06 


2 


.73 


X 


10" 


-1 





.5.5+!! 

— 


14 
20 


1 


.64 


X 


10" 


-1 


1 


.00+"'?" 

— 0.19 





.06 


8. 


.67 


X 


10" 


-1 





55+"- 

J — 


14 
20 


5 


.49 


X 


10" 


-1 


1 


00+"?" 

— 0.19 





.05 


2 


.53 


x 


10" 


-4 





00+" 


50 
()0 


1 


.58 


X 


10" 


-4 





00+"'^" 

W — 0.00 





.05 


7. 


.95 


X 


10" 


-4 


1 


.00+" 

— 


00 
50 


5 


.09 


X 


10" 


-4 





■00+"'f" 

— 0.00 





.05 


2 


.40 


x 


10" 


-3 


1 


00+" 


00 
50 


1 


.57 


X 


10" 


-3 





00+"'^ 
W — 0.00 





.05 


7. 


.89 


X 


10" 


-3 


1 


00+" 


00 
50 


5 


.03 


X 


10" 


-3 





co+0.08 
°° — 0.08 





.05 


2 


.50 


X 


10" 


-2 


1 


.00+" 

— 


00 

50 


1 


.60 


X 


10" 


-2 


1 


■00+"' 1 "! 
w — 0.50 





.05 


7. 


.84 


X 


10" 


-2 


1 


.00+" 

— 


00 
50 


1 


.93 


X 


10" 


-2 


1 


.00+"''"! 

— 0.50 





.05 


2 


.56 


X 


10" 


-1 


1 


.00+" 

— 


00 
50 


1 


.57 


X 


10" 


-1 


1 


■00+"' 1 "! 
w — 0.50 





.05 


8. 


.01 


X 


10" 


-1 


1 


.00+" 

— 


00 
50 


1 


.87 


X 


10" 


-1 


1 


.00+!!''"! 

w — 0.50 





.06 


3 


.43 


X 


10" 


-5 





.00+" 

— 


00 

00 


2 


.25 


X 


10" 


-5 





.00+"!!!; 

— 0.00 





.15 


1 


.13 


X 


10" 


-4 





.11+" 

— 


0!) 

00 


6 


.91 


X 


10" 


-5 





no+0.07 
— 0.03 





.24 


3 


.47 


X 


10" 


-4 





.51+" 

— 


17 
22 


2 


.27 


X 


10" 


-4 





00+"', 1 ,;, 

— 0.00 





.06 


3 


.5.5 


X 


10" 


-3 


1 


.00+" 

— 


00 
0.3 


2 


.21 


X 


10" 


-3 





68+"'??, 

— 0.10 





.15 


1 


.10 


X 


10" 


-2 





.9.3+" 

— 


04 
.07 


6 


.65 


X 


10" 


-3 





7 o+0.09 
' v> — 0.09 





.24 


3 


.63 


X 


10" 


-2 





.70+" 

— 


11 

2;-i 


2 


.16 


X 


10" 


-2 





.80+"'.^ 

— 0.21 





.06 


2 


.31 


X 


10" 


-5 


1 


.00+" 

— 


00 
1 ■ 


1 


.42 


X 


10" 


-5 





.00+"'"!! 

— 0.00 





.06 


7. 


.19 


X 


10" 


-5 


1 


00+" 

— 


00 

00 


4 


.56 


X 


10" 


-5 





.00+"'"!! 

— 0.00 





.06 


2 


.33 


x 


10" 


-4 


1 


00+" 


00 
00 


1 


.43 


X 


10" 


-4 





04+"'"* 
' u *-0.03 





.06 


7. 


.17 


x 


10" 


-4 


1 


00+" 

'""-0 


00 
00 


4 


.45 


X 


10" 


-4 





71+"'°? 

' ' -0.13 





.06 


2 


.23 


x 


10" 


-3 


1 


00+" 


00 
00 


1 


.45 


X 


10" 


-3 





81+"'"? 

o± -0.11 





.06 


7. 


.10 


x 


10" 


-3 


1 


00+" 


00 
00 


4 


.32 


X 


10" 


-3 





81+"'°? 

o± -0.11 





.06 


2 


.35 


x 


10" 


-2 


1 


00+" 

UU_Q 


00 
00 


1 


.43 


X 


10- 


-2 





0-1 +0.09 
' ol -0.11 





.06 


7. 


.38 


x 


10" 


-2 


1 


00+" 


00 
00 


1 


.51 


X 


10" 


-2 





81+"'°? 

o± -0.11 





.07 


2 


.53 


x 


10" 


-6 





00+" 

UU_Q 


11 

00 




0.00 x 10° 


_1 


QQ+1.00 

uu -1.00 





.07 


7. 


.55 


x 


10" 


-6 





00+" 


11 

00 


4 


.66 


X 


10- 


-6 





00+"'" 

uu -0.00 





.07 


2 


.37 


x 


10" 


-5 





10+" 


17 
05 


1 


.43 


X 


10- 


-5 





00+"" 
■ u -0.00 





.07 


7. 


.69 


x 


10" 


-5 


1 


00+" 
""-0 


00 
00 


4 


.63 


X 


10" 


-5 





00 +t).ii 

uu -0.00 





.07 


2 


.41 


x 


10" 


-4 


1 


00+" 


00 

' ' 


1 


.52 


X 


10" 


-4 





O0+0.11 
uu -0.00 





.07 


7. 


.67 


x 


10" 


-4 


1 


00+!! 


00 
00 


4 


.54 


X 


10- 


-4 





00+!!-, 1 ,, 1 , 

uu -0.00 





.07 


2 


.46 


X 


10" 


-3 


1 


.00+° 


00 
00 


1 


.43 


X 


10- 


-3 





.oo+s-iS 

— u.uu 





.07 


7. 


.48 


X 


10" 


-3 


1 


.00±° 


00 
00 


4 


.42 


X 


10- 


-3 





00+0.11 

uu -0.00 





.07 


2 


.24 


X 


10" 


-2 


1 


.00+", 


00 
00 


1 


.48 


X 


10- 


-2 





00+0.11 

uu -0.00 





.06 


2 


.43 


X 


10" 


-6 





.00+.°; 


07 
00 




0.00 x 10° 


-1 


00+100 

uu -1.00 





.06 


7 


.10 


X 


10" 


-6 





.00+", 


07 
00 


4 


.48 


X 


10- 


-6 





00+0.07 

uu -0.00 





.06 


2 


.16 


X 


10" 


-5 





.89+.°; 


OK 

11 


1 


.35 


X 


10" 


-5 





00+0.07 

uu -0.00 





.06 


7 


.16 


X 


10" 


-5 


1 


.00+°, 


Oil 

00 


4 


.33 


X 


10- 


-5 





04+0.08 
U4 -0.03 





.06 


2 


.16 


X 


10" 


-4 


1 


•ool° ; 


00 
00 


1 


.37 


X 


10- 


-4 





56+"' 11 
oo -0.11 





.06 


7 


.05 


X 


10" 


-4 


1 


.00+°, 


00 
00 


4 


.43 


X 


10- 


-4 





56+"' 11 
oo -0.11 





.06 


2 


.20 


X 


10" 


-3 


1 


•ool° ; 


00 
00 


1 


.40 


X 


10- 


-3 





56+"' 11 
oo -0.11 





.06 


6 


.79 


X 


10" 


-3 


1 


.00+° 


00 
00 


4 


.28 


X 


10- 


-3 





56+"' 11 
oo -0.11 





.06 


2 


.14 


X 


10" 


-2 


1 


•ool° ; 


00 
00 


1 


.36 


X 


10- 


-2 





56+"' 11 
OD -0.11 



Note. — Column (1) gives the cluster name. Column (2) gives the median radius, 
scaled to the virial radius of the cluster, of galaxies that go into the bin. Columns 
(3) and (5) give the median observed frame fluxes in the 8fim and 24/im channels, 
respectively, of the model SEDs that make up each bin. Fluxes are ca l culat ed by 
integrating model SEDs with random combinations of the [Asscf ct al. (201Qi) star 
forming templates across the published instrument response functions. If the SFR 
inferred from the rest frame luminosities in the model SEDs arc outside the range 
10~ 2 < SFR/(1 Mq yr" 1 ) < 10 2 the associated fluxes arc not included in the sam- 
ple. Because the Asscf ct al. (2010) templates arc not constructed to have identical 
SFRs in the S/xm and 24/xm channels, this sometimes means that an SED with a 
valid SFR in one channel will not appear in another. When a flux bin is occupied 
in one channel and not in another, the empty channel has f u — and C\ = — 1. 
This is the case for the first flux bin in AC114. Columns (4) and (6) give the MIR 
completeness (C ' \) as defined in ^4.2| The complete table is available from the elec- 
tronic edition of the journal. A sample is shown here for guidance regarding form 
and content. 
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TABLE 3 

Cluster Member Summary 



Name 


RA 


Dec 


z 


M* [1O 1O M ] 


SFR [M Q yr- 1 ] 


8 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


a3128-001 


03:30:37.7 


-52:32:57.7 


0.063 


3.2 ± 0.7 ± 1.9 


< 0.33 


0.89 


a3125-001 


03:27:20.2 


-53:28:34.6 


0.062 


3.3 ±0.7 ±2.0 


5.1 ±0.7 


3.04 


a644-005 


08:17:25.8 


-07:33:42.5 


0.071 


0.1 ±0.0 ±0.1 




1.07 


a2104-001 


15:40:07.6 


-03:17:06.8 


0.153 


3.1 ± 1.0 ± 1.9 


< 0.20 


1.33 


al689-004 


13:11:29.5 


-01:20:27.7 


0.183 


67.1 ±18.2 ±22.9 


< 1.22 


1.47 


a2163-001 


16:15:25.8 


-06:09:26.4 


0.200 


< 1.4 


< 0.50 


1.99 


msl008-001 


10:10:34.1 


-12:39:52.7 


0.308 


5.3±1.7±1.6 


< 0.87 


0.81 


acll4-001 


22:58:52.3 


-34:46:47.9 


0.303 


< 0.1 


< 0.39 


1.20 



NOTE. — T he properties of cluster member galaxies, determined using the methods described 
bv lAtle e et al. | 1120111) ■ (1) The name of this object, which is identical to the name given in Table 
2 of ^tlee^tHilT 11201 If ). (2-3) Positions of each object in J2000 coordinates, as determined from 
the ii-band images of identified cluster members. (4) Rcdshifts of each object, as determined 
by Martini ct al. (2006,2007) where available, or from the literature otherwise. (5) Stellar masses 
derived using mass-to-light ratios appropriate for each gala xy's color and assumi ng a scaled Salpctcr 
IMF with Bruzual & Chariot population synthesis model (Bell & dc Jons 2001, Table 4). The first 
uncertainty quoted gives the statistical error, and the second gives the systematic error. Systematic 
errors arc derived by applying the M/L coefficients appropriate for a Salpctcr IMF and the Pegase 
population synthesis model. (6) Star- format ion rates either from the 8/i.m luminosity, the 24/im 
luminosity or by taking the geometric mean of the two, depending on the measurements available. 
Uncertainties include only statistical errors, and upper limits ar e quoted at 3<r in the more se nsitive 
of the 8/im and 24/i.m bands. (7) Substructure parameter of Dressier &: Shcctman (1988), <5 — 
(11/<t 2 ) X [(vi oca i — v) 2 + (c i oca i — ct) 2 ] , where the local average velocities (v) and velocity dispersions 
(cr) are calculated over the 10 nearest neighbors of each galaxy. The complete table is available from 
the electronic edition of the journal. A brief sample is shown here for guidance regarding form and 
content. 
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TABLE 4 
Partial Correlation Results 







Partial r s 


Prob. 




(1) 


(2) 


(3) 


(4) 




SFR 


M» 


+0.091 


5.15 x 10" 


-01 


SFR 


-R/ R'2oo 


+0.342 


2.11 x 10" 


-05 


SFR 


s 


-0.101 


5.26 x 10" 


-01 


SFR 


s 


+0.018 


8.40 x 10" 


-01 


M» 


R/R200 


-0.087 


5.14 x lO - 


-01 


M, 


s 


-0.024 


7.93 x 10" 


-01 


M* 


s 


-0.012 


8.90 x 10" 


-01 


R/R200 




+0.188 


2.77 x 10" 


-02 


R/R200 


s 


-0.576 


1.96 x 10" 


-17 


6 


s 


+0.068 


5.41 x 10" 


-01 


Note. 


— Partial 


correlation 


results for star 



forming galaxies derived from the Spearman cor- 
relation coefficients for the variables listed in 
columns (1) and (2). Column (3) gives the strength 
of the correlation between the two variables with 
the other parameters held fixed. Column (4) gives 
the probability that a correlation at least as strong 
as that observed might occur by chance among in- 
trinsically uncorrelated data. 



